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ABSTRACT 

We present deep long-slit optical spectra for a sample of 36 Ultraluminous Infrared 
Galaxies (ULIRGs), taken with the William Herschel Telescope (WHT) on La Palma 
with the aim of investigating the star formation histories and testing evolutionary 
scenarios for such objects. Here we present the sample, the analysis techniques and a 
general overview of the properties of the stellar populations; a more detailed discus- 
sion will be presented in a forthcoming paper. Spectral synthesis modelling has been 
used in order to estimate the ages of the stellar populations found in the diffuse light 
sampled by the spectra in both the nuclear and extended regions of the target galax- 
ies. We find that adequate fits can be obtained using combinations of young stellar 
populations (YSPs, tysp l£ 2 Gyr), with ages divided into two groups: very young 
stellar populations (VYSPs, Wysp i£ 100 Myr) and intermediate-young stellar pop- 
ulations (lYSPs, 0.1 < tiYSP < 2 Gyr). Our results show that YSPs are present at 
all locations of the galaxies covered by our slit positions, with the exception of the 
northern nuclear region of the ULIRG IRAS 23327-^2913. Furthermore, VYSPs are 
presents in at least 85% of the 133 extraction apertures used for this study, being 
more significant in the nuclear regions of the galaxies. Old stellar populations (OSPs, 
^OSP > 2 Gyr) do not make a major contribution to the optical light in the majority 
of the apertures extracted. In fact they are essential for fitting the spectra in only 5% 
(7) of the extracted apertures. The estimated total masses for the YSPs (VYSPs -I- 
lYSPs) are in the range 0.18x10^° < Mysp < 50 xW^° M©. We have also estimated 
the bolometric luminosities associated with the stellar populations detected at optical 
wavelengths, finding that they fall in the range 0.07 x 10^^ < Lboi < 2.2 x 10^^ Lq. 
In addition, we find that reddening is significant at all locations in the galaxies. This 
result emphasizes the importance of accounting for reddening effects when modelling 
the stellar populations of star- forming galaxies. 
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1 INTRODUCTION 

The launch of the Infrared Astronomical Satellite (IRAS) in 
1983 triggered a development of infrared astronomy that is 
still ongoing. Although the majority of the sources discov- 
ered by IRAS were modest infrared emitters, some objects 
had mid- to far-infrared (MFIR) luminosities comparable 
to the bolom e tric luminos it ies of o ptically -se lected QSOs 
I Houck et al.1 1 1984 Il985l : ISoifed . Il984al lbl: ISoifer et al.l . 



1 19871 '). Such galaxies are classified as Luminous (Li,. > 10^ 



Lq) or Ultraluminous {Li,. > 10^^ L0) infrared galaxies 
(LlRGs/ULlRGs) depending on their MFIR luminosities. 

The prodigious infrared emission of ULIRGs is gener- 
ally attributed to the optical/UV light of luminous central 
sources reprocessed by dust. Several infrared spectroscopic 
studies have provided evidence that the power source of 
ULIRGs has a compo site nature, with a mixture of AGN 
and starburst activit y ([Genzel et al.','l998VLutz et al.','l999|; 
llm; 
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Armus et al. I, I2OO7I : llmanishi et a l.. 2007 : Farrah ct ahl. 
2007^ . Similar conclusions are reache d in the relativ e ly few 
studi e s carried out at U V/optical l|Veilleux et al.l, 19951 , 
19991 : iFarrah et al.1 . 120031 ) and X-ray (|Franceschini et al.l . 
200ll . l2003l ') wavelengths. Therefore these objects provide us 
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with an excellent oportunity to study the links between AGN 
and starburst phenomena. In addition, ULIRGs are almost 
unv ariably as sociated with galaxy mergers and interactions 
fe.g. lKim efal .. 2002. : . Veilleux ct al. , ,2002i ). Thus, they also 
represent ideal objects to test models of galaxy evolution via 
major galaxy mergers (e.g. Mihos & Hernquist, 1996; Barnes 
& Hernquist, 1996; Springel et al., 2005). 

Stellar population studies have the potential to provide 
key information about the histories of the merger events 
asscoaited with ULIRGs. To date, two complementary 
approaches have been used to investigate the stellar popu- 
lations in such objects: one based on photometric analysis 
of images of the galaxies taken at different wavelengths, 
focussed mainly on studying the young stellar popu- 
lation s associated with t he brig ht knots (|Surace et al.l . 



19981: [S r 



& Sanders 
Surace fc Sanders. 2000; 



s png ri 
1 19991: 



20001 : 



ISurace e t al] , 
Wilson et a l.. 2006J; the other 
comprising spectroco p ic ana lys is of UV/optical spectra 



comprising spectroco p ic ana lys is or u v / optical spectra 
llCanalizo fc Stocktonl |2000allbl l200ll: iFarrah et all. 12005 



2007i : iRodriguez Zaurm et al. 
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20081 ). Overall, there have been relatively few studies of 



the stellar populations in ULIRGs. The existing studies 
mentioned above either concern relatively small samples of 
ULIRGs, or concentrate on high surface brightness regions 
which may give a misleading impression of the stellar 
populations of the systems as a whole. 

To remedy this situation we have obtained deep, wide 
spectral coverage, intermediate-resolution long-slit spectro- 
scopic observations for a sample of 36 ULIRGs with redshifts 
z < 0.18 using the ISIS spectrograph on the WHT. With the 
aim of studing in detail the properties of the stellar popu- 
lations in the ULIRGs in our sample, and thereby testing 
evolutionary scenarios, we model the spectra extracted from 
a series of apertures using a large number of different com- 
binations of stellar popula tions synthesis model templates 
IIBruzual fc Charlotl . |2003| ). In this paper we present the 
sample and the spectral synthesis modelling results, while 
a more detailed analysis of the results is presented in the 
forthcoming paper (Rodriguez Zaurm et al., in prep). De- 
tailed analyses of two of the objects in the sample, IRAS 
13451-^1232 (PKS1345+12) and IRAS 15327-^2340 (Arp 
220) were presented in Rodriguez Zaurm et al. (2007, here- 
after RZ07) and Rodriguez Zaurm et al. (2008, hereafter 
RZ08) respectively. 

Throughout this paper, we assume a cosmology with Ho 
= 71 km s"^ Mpc"\ Q.0 = 0.27, D.a = 0.73. AU the values 
presented in this paper have been adapted to this cosmology. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The Sample 

In the first place, we observed a complete RA- and 
declination-limited sub-sample of the lKim fc Sanders! (|l998l ) 
1 Jy sample of ULIRGsQ, with RAs in the range 12 < RA < 
1:30 h, declinations S > -23 degrees, and redshifts z < 0.13. 



The redshift limit was chosen to ensure that the objects are 
sufficiently bright for spectroscopic study, and also to keep 
the size of the sample tractable for deep observations on 4m 
class telescopes. Throughout the paper, we will refer to this 
sample of 26 objects as the complete sample (CS). Table [1] 
summarizes some of the properties of the ULIRGs in the CS. 

As shown in the table, only 8 ULIRGs with warm mid- 
to far-IR colours (f25/f60 > 0.20 are included in the CS. 
However, one of the aims of the project was to study the 
presence of an evolutionary link between cool and warm 
ULIRGs. Therefore, we extended the sample to include 6 
higher-redshift objects classified as warm ULIRGs: IRAS 
13305-1739, IRAS 14252-1739, IRAS 16156+0146, IRAS 
17044+6720, IRAS 23060+0505 and IRAS 23389+0303. 
In addition, since we also aim to test the evolutionary 
scenario in which ULIR Gs evolve into QSOs or a radio 
galaxies (e.g. Surace et al. . 1998; Surace & Sa ndersl Il999l : 
ICanalizo fc Stocktonll200ll : [Tadhunter et al.1 . 120051 '). we de- 



cided to include IRAS 17179+5444, which is classified as a 
cool ULIRG with a Sy2-like optical spectrum. Finally, three 
ULIRGs were included in the sample (IRAS 08572+3915, 
IRAS 10190+1322 and IRAS 10494+4424) which have red- 
shifts z < 0.13, but fall outside the RA range of the CS. One 
of them, IRAS 08572+3915, is classified as a warm ULIRG. 
Thus, the total sample discussed in this paper comprises 36 
objects: 21 cool ULIRGs, and 15 warm ULIRGs. We will 
refer to this sample as the extended sample (ES). With an 
upper limiting redshift of z < 0.18, the ES includes all the 
warm ULIRGs within the RA and declination range of the 
CS with redshifts z < 0.15. Table [2] summarizes some of 
the properties of the 10 ULIRGs which, together with the 
CS, comprise the ES. Optical and near-IR images fo r all the 
objects in the ES are presented in lKim et al.l l|2002h . 

Figure [1] presents histograms comparing the full IJy 
sample, the CS and the ES in terms of IR-luminosity, f25/f60 
flux ratio and morphological type distributions, respectively. 
The upper panel shows that both the CS and the ES are rep- 
resentative of the IJy sample in terms of IR-luminosities, 
while the middle panel of the figure shows that, as ex- 
pected, the ES is biased towards objects classified as warm 
ULIRGs (f25/f60 > 0.2). Using the morphologi cal classifi- 
cation scheme described in lVeilleux et al.l (|2002l '). the lower 
panel in Figure [1] shows the distribution of the different mor- 
phological type s of the o bjects comprising each sample. To 
summarize (see lVeilleux et al. . 2002. for further details), the 
different classes are: cIeiss I (first approach), class II (first 
contact), class III (pre-merger), class IV (merger), class V 
(old merger), Tpl (triple) and Iso (Isolated). As is clear in 
the figure, both the CS and the ES samples account for 
most of the morphological types of objects found in the full 
1 Jy sample, with the exception of those classified as 
This is due to the fact that all but 4 of the IVa sources in 
the 1 Jy sample are at redshifts z > 0.13 (one of these 4 
sources, IRAS 14060+2919, is included in the CS), and the 
objects in the CS are selected at low (z < 0.13) redshifts. 



^ This is in itself a complete flux-limited sample of 118 ULIRGs 
identifled from the IRAS Faint Source Catalogue (FSC) and se- 
lected to have a 60/im flux densities greater than 1 Jy in a region 
of the sky with 5 > -40° and |fe| > 30°. 



^ The quantities f25 and /go represent the IRAS flux densities in 
Janskys at 25Mm and GO^tm. 

^ In the IVeilleux et al. morphological classification, ob- 

jects classified as IVa are those with prominent tidal features and 
diffuse extended central regions, but only one nucleus. 
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Object Name 


z 


RA 


DEC 


log 


£25 / fgo 


Nuclear 


NS 


IC 


ST 


IRAS 




(J2000.0) 


(J2000.0) 






structure 


(kpc) 














(Lq) 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


00091-0738 


0.118 


00 11 43.4 


-07 22 06 


12.23 


0.08 


Double 


2.2 


Illb 


H 


00188-0856 


0.128 


00 21 26.0 


-08 39 29 


12.37 


0.14 


Single 


- 


V 


L 


01004-2237 


0.129 


01 02 51.2 


-22 21 51 


12.28 


0.28 


Single 


- 


V 


H 


12072-0444" 


0.129 


12 09 45.4 


-05 01 14 


12.39 


0.21 


Double 


2.8 


IVb 


S2 


12112+0305 


0.073 


12 13 47.3 


02 48 34 


12.32 


0.06 


Double 


4.0 


Illb 


L 


12540+5708 


0.042 


12 56 15.0 


56 52 17 


12.54 


0.27 


Single 


- 


IVb 


SI 


13428+5608" 


0.037 


13 44 41.8 


55 53 14 


12.14 


0.10 


Double 


0.7 


IVb 


S2 


13451+1232 


0.122 


13 47 33.3 


12 17 24 


12.32 


0.35 


Double 


4.2 


Illb 


82 


13539+2920 


0.108 


13 56 10.9 


29 05 29 


12.04 


0.06 


Double 


7.4 


Illb 


H 


14060+2919 


0.117 


14 .08 17.5 


29 04 57 


12.07 


0.08 


Single 


- 


IVa 


H 


14348-1447 


0.083 


14 37 37.3 


-15 00 20 


12.32 


0.07 


Double 


5.1 


Illb 


L 


14394+5332'' 


0.105 


14 41 04.3 


53 20 08 


12.08 


0.17 


Multiple 


54.0 


Tpl 


S2 


15130-1958 


0.109 


15 15 55.6 


-20 09 18 


12.09 


0.20 


Single 


- 


IVb 


S2 


15206+3342 


0.125 


15 22 38.0 


33 31 36 


12.22 


0.20 


Single 


- 


IVb 


H 


15327+2340 


0.018 


15 34 57.1 


23 30 10 


12.21 


0.07 


Double 


0.4 


Illb 


L 


15462-0450 


0.100 


15 48 56.6 


-04 59 36 


12.21 


0.15 


Single 




IVb 


SI 


16474+3430 


0.111 


16 49 14.7 


34 25 13 


12.15 


0.09 


Double 


6.9 


Illb 


H 


16487+5447 


0.104 


16 49 47.8 


54 42 34 


12.16 


0.07 


Double 


5.7 


Illb 


L 


17028+5817 


0.106 


17 03 41.8 


58 13 48 


12.14 


0.04 


Double 


24.6 


Ilia 


L 


20414-1651 


0.086 


20 44 17.4 


-16 40 14 


12.18 


0.08 


Single 




IVb 


H 


21208-0519 


0.130 


21 23 28.7 


-05 06 59 


12.05 


0.13 


Double 


14.9 


Ilia 


H 


21219-1757 


0.112 


21 24 42.5 


-17 44 40 


12.10 


0.42 


Single 




IVa 


SI 


22491-1808 


0.076 


22 51 49.0 


-17 52 27 


12.13 


0.10 


Double 


2.4 


Illb 


H 


23233+2817 


0.114 


23 25 48.7 


28 34 19 


12.04 


0.22 


Single 




Iso. 


S2 


23234+0946 


0.128 


23 25 56.0 


10 02 52 


12.09 


0.05 


Double 


7.8 


Illb 


L 


23327+2913 


0.107 


23 35 12.5 


29 30 05 


12.10 


0.10 


Double 


24.0 


Ilia 


L 


Table 1. The Complete sample 


of 26 ULIRGs with redshifts z < 0.13 and declinations S > -23 


degrees. Col (1) 


: object 


the IRAS Bright Galaxy survey dSoifer 


et al.Ul986l,ll9S 


TI). Col (2): optical redshifts from|Kim & 


SandersI 1 


1998|). 


Cols (3) 



ascension (hours, minutes and seconds) and declination (degre es, arcminutes and arc seconds) of the IRAS source position as listed in the 
Faint Source Database (FSDB). Col (5): IR luminosity from lRim &: SandersI l ll998l ). Col (6): f25/f60 mid- to far-infrared colour ratio, 
where f25 and fgo represen t the I RAS flux densities (non-colour c orrected) in units of Jy at 25 and 60 fj,m respectiyely. C ols (5) and (6): 
both f rom lKim &: Sander^ |l99^. Col (7): nuclear structure (from llCim et alll2002l) . Col (8): nuclear separation (NS) from lVeilleux et al.l 
||2002| ). Col (9): interaction class. Cl ass I: first approach; Class II: first contact; Class III: pre-merger; Class IV: merger; Class V: old 
merger; Tpl: Triple; Iso: isolated, (see Veilleux et al.L [20021 . for details). Col (10): optical nuclear spectral type (from fVeilleux et al.L[l995l : 
iKim fc Sander£l.1l99^;IVeilleux et allll999h . H: HII; L: LINER; SI: Seyfert 1; S2: Seyfert 2. 

Kim et al. I l l2002li" classi fied these objects as single nucleus systems. Howeyer, a double nucleus structure is reyealed in the studies of 



iDasyra et ahl ll2006a^ and lScoyille et al.l 1 I2OOOI) for IRAS 12072-0444 and IRAS 13428+5608 respectiyely. These objects will be classified 
as double nucleus systems for the work presented here. 

^ This source is a multiple (> 2 nucleus) system. The two main components are separated 54 kpc, but the eastern component itself 
comprises two close nuclei. 
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Nuclear 
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ST 
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(J2000.0) 


(J2000.0) 






structure 


(kpc) 














(L©) 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


08572+3915 


0.058 


09 00 25.0 


39 03 56 


12.15 


0.22 


Double 


6.0 


Illb 


L 


10190+1322 


0.077 


10 21 41.9 


13 07 01 


12.04 


0.11 


Double 


5.9 


Illb 


H 


10494+4424 


0.092 


10 52 22.2 


44 08 59 


12.17 


0.04 


Single 




IVb 


L 


13305-1739 


0.148 


13 33 15.2 


-17 55 01 


12.25 


0.33 


Single 




V 


S2 


14252-1550 


0.149 


14 28 01.4 


-16 03 43 


12.19 


0.20 


Double 


8.8 


Illb 


L 


16156+0146 


0.132 


16 18 08.2 


01 39 21 


12.08 


0.24 


Double 


8.0 


Illb 


S2 


17044+6720 


0.135 


17 04 28.5 


67 16 34 


12.17 


0.28 


Single 




IVb 


L 


17179+5444 


0.147 


17 18 55.1 


54 41 50 


12.24 


0.14 


Single 




IVb 


S2 


23060+0505 


0.173 


23 08 34.2 


05 21 29 


12.48 


0.37 


Single 




IVb 


S2 


23389+0303 


0.145 


23 41 31.1 


03 17 31 


12.13 


0.28 


Double 


5.2 


Illb 


S2 



Table 2. Same as Table 1, but for the 10 additional objects included in the extended sample (ES) discussed in this paper. 
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Figure 1. Upper panel: Histograms sliowing the IR luminosity distributions for the full 1 Jy sample of iKim fc Sander^ lll998l ). the 
complete sample (CS) and the extended sample (ES). The figure shows that both the CS and the ES are representative of the IJy 
sample in terms of IR luminosities. Middle panel: same as the upper panel, but for the f25/f60A'™ ratio, used to classifiy the ULIRGs 
as cold (f25/f60 < 0.2) or warm (f25/f60 ^ 0.2) objects. The figure shows that the CS is representative of the IJy sample in terms of 
the fgs/feo-ratio distribution. On the other hand, as expected, the ES is biased towards warm objects. Lower panel: same as the upper 
panel, but for the interaction class distribution, using the morphological classification scheme described in Veilleux et al. (2002). The 
figure shows that the CS and the ES both have similar distributions of objects in interaction class to the IJy sample, with the exception 
of the IVa class (see the text for an explanation). 



In addition, among the 23 objects classified as IVa in the 
IJy sample, only one is classified as a warm ULIRG, and 
therefore, these objects were not selected, a priori, to be ob- 
served and included in the ES. To summarize, the CS and 
the ES samples are representative of ULIRGs in the local 
universe, although the ES contains a higher proportion of 
warm objects. 

2.2 Optical spectroscopy 

Long-slit spectra of all the objects in the ES were taken 
in July 2005, December 2005 and May/July 2006, with the 
ISIS dual-beam spectrograph on the 4.2-m William Herschel 



Telescope (WHT), on La Palma, Spain. With the exception 
of IRAS 22491-1808, described below, we used the R300B 
grating with the EEV12 CCD, and the R316R grating with 
the MARC0NI2 CCD on the blue and the red arm respec- 
tively, along with a dichroic cutting at ~ 5300 A. A 2x2 
binning mode was used during the observations, leading to 
a final spatial scale of 0.4 arcsec/pix for both arms, and a 
dispersion of 1.72 A/pix in the blue and 1.65 A/pix in the 
red arm. The observed useful wavelength range is ^ 3300 - 
7800 A. A slit of width 1.5 arcsec was used for all the objects 
in the ES. In the case of IRAS 22491-1808, the settings in 
the blue arm were the same as described above. However, on 
the red arm, the R600R grating with the MARC0NI2 CCD 
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was used, with 2x2 binning mode. The final spatial scale 
was 0.4 arcsec/pixel and the dispersion was 0.89 A/pixel. 
The observed useful wavelength range is ~ 3300 - 7100 A. 
Details of the observations for all the objects included in the 
ES are presented in Table [3] 

The position angles (PA) of the slit were carefully se- 
lected to cover the main features of the objects (nuclei,tidal 
tail, bridges, knots). In order to minimize the effects of dif- 
ferential atmospheric refraction, the objects were observed 
either along parallactic angle or at small airmass (AM < 
1.1). For those objects with a double nucleus structure, 
the slits covered both nuclei, with the exceptions of IRAS 
08572+3915 and PKS 1345-f 12, for which it was not possi- 
ble to observe the object at the parallactic angle and we did 
not have sufficient time to observe both nuclei during the 
night of the observations. In the cases of the double nucleus 
objects IRAS 16487-^5447 and IRAS 17028-^5817, because 
of the parallatic angle constraints, the two nuclei in each 
system were observed with separate slit PAs. 

The data were reduced (bias subtracted, flat fleld cor- 
rected, cleaned of cosmic rays, wavelength calibrated and 
flux calibrated) and straightened before extraction of the in- 
dividual spectra using the standard packages in IRAF and 
the STARLINK packages FIGARO and DIPSO. The wave- 
length calibration accuracy, measur ed as the mean shift be - 
tween the measured and published (|Qsterbrock et all . Il996l ) 
wavelength of night-sky emission lines, is ~ 0.35 A for the 
blue spectra and ~ 0.25 for the red spectra. No significant 
variations of these values were found for the observations of 
the objects comprising the sample used in this paper. The 
spectral resolutions, calculated using the widths of the night- 
sky emission lines (FWHM), are in the ranges of 3.4 - 6.03 
A and 4.44 - 5.84 A for the blue and red arms, respectively. 

Accurate relative flux calibration is vital for the pur- 
pose of this project, and therefore several (4 - 6) photo- 
metric standard stars were observed during each night of 
the observations. The relative flux calibration accuracy is 
estimated to be ± 5% over the entire useful spectral range 
(^3300 - 7800 A). This is confirmed by the excellent match 
between the fluxes in the overlapping region of the blue and 
red arm spectra. 



3 ANALYSIS 

3.1 Aperture extraction 

The relatively low redshifts of the objects in the ES, the wide 
spectral coverage, intermediate spectral resolution, excellent 
flux calibration and high S /N of the spectra, make it possi- 
ble to perform detailed studies of their stellar populations. 
The extraction apertures were selected from spatial cuts of 
the 2-D frames in the line-free continuum wavelength range 
4400 - 4600 A, based on the visible extended structures and 
the requirement that the apertures are large enough to have 
a sufficiently high S/N ratio for further analysis. In order 
to compare the stellar populations between the objects in 
our sample, apertures with a metrical scale of 5kpc cen- 
tred on the main nuclei were extracted for all the objects in 
the extended sample (ES), including separate extractions for 
multiple nuclei in individual sources. As mentioned before, a 
detailed analysis of the particular cases of IRAS 1341-M232 



(PKS 1345+12) and IRAS 15327+2740 (Arp 220) can be 
found in RZ07 and RZ08 respecively. In the case of PKS 
1345+12 the aperture labelled as 'NUC AP' in the RZ07 
paper has a metric scale of 5.2 kpc and, hereafter, we will 
refer to this as the 5kpc aperture for this galaxy. In the case 
of Arp 220, the closest object in the sample (z = 0.018, re- 
sulting in a scale of 0.363 kpc arcsec"'^ and a distance of 
77.6 Mpc) , a 5 kpc aperture covers a large region of the slit 
and, therefore, was not selected a priori when performing the 
study presented in RZ08. In order to compare Arp 220 with 
the other ULIRGs included in the sample used for the work 
presented here, we used the slit with PA 160 and extracted 
a 5kpc aperture, as shown in Figure (2] sampling the central 
region of the galaxy (see Figure 1 in RZ08). A second set of 
apertures was then selected to sample the spatial features 
of those objects in the ES showing tails, bridges and other 
diffuse structures. A total of 133 apertures was extracted for 
the ULIRGs in the ES. Figure [5] shows the spatial profiles of 
the 2-D frames in the wavelength range 4400 - 4600 A along 
with the extraction apertures. 

The extracted spectra, in the rest frame, are shown in 
Figure |6l In the particular case of Arp220 we extracted a 
total of 24 apertures. As discussed in RZ08, the morpholo- 
gies of the extracted spectra (and the properties of the stel- 
lar populations) are similar across the entire extent of the 
galaxy. Therefore, in Figure (6] we decided to show a rep- 
resentative sample of the whole set of extracted apertures. 
These are the apertures labelled as APtotal in Figure [2] 
for each of the slit PAs used for this galaxy (see RZ08 for 
more details). 

3.2 Modelling technique 

To study in detail the properties of the stellar populations 
we have modelled the full spectral range of the extracted 
spectra using the stellar populati on synthesis model tem- 
plates of lBruzual fc Chariot! (120031 ) . These spectra were cre- 
ated assuming an instantaneo us burst o f star f ormation, so- 
lar metallicity and using the ISalpeteil (|l955l ) initial mass 
function (IMF) with lower and upper mass cut offs: m_L = 0.1 
Mq and mu = 100 Mq (see lTadhunter et al.1 . [20051 . for a jus- 
tification of these assumption) . On the other hand, ULIRGs 
have large gas and dust contents. Therefore, it is important 
to carefully account for reddening eff ects. Reddened syn - 
thetic spectra were created using the iCalzetti et al.1 (|2000| ) 
red dening law, ap propriate for starburst galaxies, and also 
the ISeatonI (|l979l ) reddening law, representing the Galactic 
extinction case. During this process, the intrinsic extinction 
was modelled using a foreground screen geometry. Overall, 
there is good consistency between results obtained with the 
two reddening laws, suggesting that the main results for the 
ULIRGs in the CS and the ES samples are not sensitive to 
the details of the reddening law assumed at optical wave- 
lengths 

Of the 133 extraction apertures, four sample the nuclear 
regions of Syl galaxies (5kpc and Ap D for Mrk 231, 5kpc for 
IRAS 15462-0405, and 5kpc in the case of IRAS 21219-1757) 
and, therefore, it was not possible to model the extracted 
spectra due to the strong AGN contamination. In addition, 
Ap A for IRAS 23233+2817 samples a region to the south 
of the nucleus, includin g the bright kno t detected in the 
K-band image shown in Kim et all (|2002l ') (their Figure 1). 
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Object name 


PA 


Arm 


Exposure 


Airmass 


Seeing FWHM 


IRAS 







(s) 




(arcsec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




00(191-0738 





B 


2700 


1.24 


1.2 








R 


2700 


1.24 


1.2 




OOlRS-OS'ifi 

\J\J -LfJO \J(JtJ\J 


332 


B 


2100 


1.33-1.37 


1.3 








R 


2100 


1.33-1.37 


1.3 




01 004-22*^7 


352 


B 


2100 


1.60 


2.8 








R 


2100 


1.60 


2.8 






265 


B 


3600 


1.14-1.32 


1.3 








R 


3600 


1.14-1.32 


1.3 




iniQn-i-i'i22 


64 


B 


3600 


1.04 


1.2 








R 


3600 


1.04 


1.2 




104944-4424 


167 


B 


3600 


1.03-1.07 


1.2 








R 


3600 


1 03-1.07 


1.2 




12072-0444 


45 


B 


2700 


1.30-1.42 


1.1 








R 


2700 


1.30-1.42 


1.1 




1211 2-OSO^ 


15 


B 


2700 


1.11 


1.2 








R 


2700 


1.11 


1.2 






180 


B 


2700 


1.13 


0.7- 


1.6 






R 


2700 


1.13 


0.7- 


1.6 


13305-1739 


342 


B 


2700 


1.53-1.62 


1.7 








R 


2700 


1 53-1 62 


1.7 




1 3428-t-5fi08 


180 


B 


2700 


1.12 


1.7 








R 


2700 


1.12 


1.7 




1 3451 -t-1 232 


160 


B 


3600 


1.05 


1.3 








R 


3600 


1.05 


1.3 






230 


B 


3600 


1.1 


1.7 








R 


3600 


1.1 


1.7 




1 3'S3Q-4-2Q20 


285 


B 


2700 


1.01 


1.1 - 


1.6 






R 


2700 


1.01 


1.1 - 


1.6 


14060-1-2919 


83 


B 


2700 


1.06-1.11 


0.9 - 


1.3 






R 


2700 


1.06-1.11 


0.9 - 


1.3 




40 


B 


2700 


1 63-1 82 

X X,Oij 


0.8 








R 


2700 


J. .UO _L . Oil 


0.8 




14348-1447 

XrrOrrO X rrrr 1 


35 


B 


2700 


1.61-1.80 


1.4 - 


2.1 






R 


2700 


1.61-1.80 


1.4 - 


2.1 


1 4394-1- "^332 


92 


B 


3600 


1.47-1.70 


1.6 








R 


3600 


1.47-1.70 


1.6 




1 51 30-1 Q58 

XtJ -L OU Xi7(JO 


358 


B 


2700 


1.53 


0.8 - 


1.1 






R 


2700 


1.53 


0.8 - 


1.1 


15206-1-3342 


90 


B 


2700 


1.23-1.28 


0.9 - 


1.8 






R 


2700 


1 23-1 28 


0.9 - 


1.8 


15327-1-2340 


160 


B 


2700 


1.03-1.07 


1.3 - 


2.0 






R 


2700 


1.03-1.07 


1.3 - 


2.0 




75 


B 


2700 


1.14-1.19 


0.6 








R 


2700 


1.14-1.19 


0.6 






75* 


B 


3600 


1.27-1.46 


0.7 








r> 


ODUU 




0.7 




15462-0450 


10 


B 


2700 


1.2 


1.6 








R 


2700 


1.2 


1.6 




16156-1-0146 


313 


B 


2700 


1.16-1.21 


1.2 - 


1.7 






R 


2700 


1.16-1.21 


1.2 - 


1.7 


16474-1-3430 


166 


B 


2700 


1.01 


1.0 








R 


2700 


1.01 


1.0 




16487-I-5447NE 


110 


B 


2700 


1.12 


1.2 








R 


2700 


1.12 


1.2 




164874-5447NW 


165 


B 


2700 


1.26-1.32 


1.1 








R 


2700 


1.26-1.32 


1.1 





Table 3. Log of the spectroscopic observations. In the case of Arp 220 (IRAS 15327H-2340), a slit in PA 75 was also used with 10 arsec 
offset to the north, with respect to the one referred as PA75 in the table: this is labelled as PA75*. The seeing was estimated from DIMM 
seeing monitor measurements. 
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Object name 
IRAS 
(1) 


PA 

o 

(2) 


Arm 

(3) 


Exposure 
(s) 
(4) 


Airmass 
(5) 


Seeing 
(arcsec) 
(6) 


17028+5817NE 


130 


B 


2700 


1.23-1.28 


1.2 








R 


2700 


1.23-1.28 


1.2 


17028+5817NW 


110 


B 


2700 


1.31-1.39 


0.9 








R 


2700 


1.31-1.39 


0.9 


170444 


-6720 


160 


B 


2700 


1.20-1.29 


1.4 - 1.9 








R 


2700 


1.20-1.29 


1.4 - 1.9 


171794 


-5444 


122 


B 


2700 


1.22-1.28 


0.9 








R 


2700 


1.22-1.28 


0.9 


20414-1651 


5 


B 


2700 


1.42 


1.4 - 1.7 








R 


2700 


1.42 


1.4 - 1.7 


21208-0519 


15 


B 


2700 


1.20 


1.2 








R 


2700 


1.20 


1.2 


21219-1757 


333 


B 


3600 


1.54-1.67 


1.5 








R 


3600 


1.54-1.67 


1.5 


22491-1808 


338 


B 


3600 


1.54-1.58 


1.0 - 1.5 








R 


3600 


1.54-1.58 


1.0 - 1.5 


230604 


-0505 


338 


B 


3600 


1.09-1.13 


0.7 








R 


3600 


1.09-1.13 


0.7 


232334 


-2817 


168 


B 


2100 


1.00 


1.2 








R 


2100 


1.00 


1.2 


232344 


-0946 


115 


B 


2700 


1.06 


1.1 








R 


2700 


1.06 


1.1 


233274 


-2913 


175 


B 


2700 


1.00 


1.7 








R 


2700 


1.00 


1.7 


233894 


-0300 


24 


B 


3600 


1.12 


0.7 








R 


3600 


1.12 


0.7 



Table 3. Continued 



Our extracted spectrum for this aperture shows that it is 
strongly contaminated by a Galactic M star (see Figure 3). 
Therefore no attempt was made to model the spectrum for 
this aperture. 

Prior to the modelling, the spectra were corrected for 
Galacti c reddening us i ng the far-IR based maps of extinc- 
tion by ISchlegel et all l|l998h . An additional potential con- 
taminant of the stellar e mission in the optical is the nebular 
continuum emission (e.g. iDickson et all Il995l ). In order to 
account for this effect, we subtracted a nebular continuum 
for those apertures for which strong emission lines were de- 
tected, following the technique described in RZ07, i.e. we 
used Ha to generate the nebular continuum and considered 
two extreme cases: (i) maximum nebular continuum assum- 
ing no reddening of the emission-line region; (ii) zero neb- 
ular continuum, corresponding to high reddening. The final 
nebular continua were subtracted from the spectra prior to 
performing the modelling. The modelling was carried out for 
both the nebular-corrected and uncorrected spectra. 

Overall, a total of 128 extraction apertures were mod- 
elle d. To perform the fit, we have used the CONFIT code 
(see iRobinson et al.l I2OO0I and RZ08 for details) which as- 
sumes two stellar components plus a power-law in some 
cases. To summarize, the CONFIT approach consists of a 
direct fit to the overall continuum shape of the extracted 
spectra using a minimum technique. For each spectro- 
scopic aperture, the flux is measured in several wavelength 
bins (typically ~50 - 70) chosen to be as evenly distributed 
in wavelength as possible, and to avoid strong emission lines 
and atmospheric absoption features. A relative flux calibra- 
tion error of 5% was assumed during the modelling. Models 



with Xred ~ 1 deemed to provi de an acceptable fi t to th e 
overall continua (see discussion in iTadhunter et al. I. I2OO5I ). 
From these, the best fitting models were selected based on 
a visual inspection of the fits to the detailed absorption fea- 
tures that have relatively little emission line contamination, 
such as high order Balmer lines. He I lines. Call K A3934, 
G-band A4305, Mgib A5173 band. In addition, for those ob- 
jects classified as Hll-galaxies, we used the Ha equivalent 
widths to further constrain the ages of the VYSlQ- 

The model fit results are quantified in terms of the per- 
centage contribution of the different stellar components in 
a normalising bin. Due to the different spectral morpholo- 
gies of the extracted spectra, it was not possible to select 
exactly the same normalising bin for all the extracted aper- 
tures. However, the normalising bin was always selected to 
be located within the wavelength range 4400 - 4800 A, and 
usually extends ~ 100 A. 

Since CONFIT allows for a maximum of two stellar 
components plus a power-law in some cases and the stellar 
contents of the galaxy are not known a priori, we used var- 
ious combinations of spectra during the modelling process 
in order to determine the properties of the dominant stellar 
populations. We would like to emphasize that the techniques 
and the modelling philosophy used in this paper are based 
on the resuhs of RZ08 for the ULIRG Arp 220, and to a 
lesser extent on those of RZ07 for the case of PKS1345-I-12. 

* Note that since the continua of the spectra are affected by di- 
lution due to the emission of older stellar components, the H« 
equivalent width measurements provide upper limits for the age 
of the VYSP. 
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IRAS 00091-0738 IRAS 00188-0856 IRAS 01004-2237 




Arcsec S Arcsec SE NW Arcsec 



IRAS 08572+3915 



IRAS 10190+1322 



IRAS 10494+4424 






IRAS 12072-0444 IRAS 12112+0305 IRAS 12540+5708 (Mrk231) 




Arcsec NW Arcaec SE s Arcsec 



IRAS 13428+5608 (Mrk273) - extended IRAS 13451 + 1232 (PKS1345+12) PA 160 IRAS 13451 + 1232 (PKS1345+12) PA 230 




NW Arcsec SE ne Arcsec SW NE Arcsec SW 



Figure 2. Spatial profiles of the 2-D spectra for tlie wavelength range 4400 - 4600 A showing the positions of the extraction apertures. 
The symbols N, S, E and W at the positive and negative ends of the x-axis represent the directions of the slit on the sky. For example, 
in the case of IRAS 00091-0738, positive x corresponds to south, negative x corresponds to north. 



The properties of the stellar populations in ULIRGs I: sample, data and spectral synthesis modelling. 



IRAS 1434B-1447 



IRAS 14394+5332 



IRAS 15130-1958 




Arcsec 
IRAS 15206+3342 




IRAS 15327+2340 (Arp220) PA 75* 

TOTAL - 




Arcsec 
IRAS 16474+3430 
Skpcl |kpcll 





Arcsec 

IRAS 15462-0450 



Arcsec 
IRAS 16487+5447 E 

5kpc 




Arcsec 
IRAS 17028+5817 W 
Skpi 




Arcsec 
IRAS 20414-1651 





Arcsec 

IRAS 15327 + 2340 (Arp220) PA 75 




Figure 2. 
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IRAS 21219-1757 IRAS 22491-1808 IRAS 23060+0505 




N Arcsec S N Arcsec 



Figure 2. 



Throughout this paper, we define young stellar popula- 
tions (YSPs) as stellar components with ages tysp< 2 Gyr, 
and old stellar populations (OSP) as components with ages 
tosp>2 Gyr. With the aim of better describing each combi- 
nation, it is convenient to divide the YSPs into two groups: 

• very young stellar populations (VYSP): stellar compo- 
nents with ages tvYSP < 0.1 Gyr; 

• intermediate-age young stellar populations (lYSP): stel- 
lar components with ages in the range of 0.1 < tiYSP < 2 
Gyr. 

Generally, for models that include an old plus a young 
and/or an intermediate stellar component it is not possible 
to distiguish between OSP ages in the range 2 < tosp < 
12 Gyr, and therefore, we decided not to use such stellar 
populations during the modelling analysis described below. 
Instead, we fixed the age of the OSP at 12.5 Gyr. The dif- 
ferent combinations used during the modelling were: 

• Combination I: as a first approach we used a two 
component model comprising a YSP (Iysp < 2 Gyr) with 
varying reddening (0.0 < E{B — V) < 2.0, increasing in 
steps of 0.1), along with an unreddenetfl OSP of age 12.5 
Gyr. Accounting for an old component covers the case in 
which one or more of the merging galaxies is early-type or 

* Although we cannot entirely rule the idea that these stellar 
populations are significantly reddened, assuming zero reddening 
for OSP is reasonable since such populations are likely to be sit- 
uated in the extended bulge/halo in the foreground of the sys- 
tem (along the line of sight), or else completely obscured by the 
circum-nuclear dust in the background of the system 



bulge-dominated galaxy. It is worth mentioning that this 
combination is almost identical to the one used in RZ07 for 
the ULIRG PKS1345-I-12. However, the version of the fitting 
code used here represents an updated version of the one used 
in that paper. Therefore, in order to compare the results ob- 
tained or this galaxy with those of the other ULIRGs in our 
sample, we decided to model again the extracted spectra 
for PKS1345-I-12 using this combination and the updated 
version of the code. The results obtained here are consis- 
tent with those of RZ07. For consistency, the spectra of this 
source were also modelled using the other two combinations 
described below. 

• Combination II: this combination consists of three 
components: an OSP of age 12.5 Gyr and zero reddening, 
along with a YSP with variable reddening (0.0 < E{B—V) < 
2.0, increasing in steps of 0.1), and a power-law {F\ oc A") 
with a spectral index in the range — 15 < a < 15. The power- 
law is included to repre sent either scattered or direct AGN 
continuum component (|Tadhunter et all 120021 ') . or a highly 
reddened VYSP. 

• Combination III: for ULIRGs as a class, there is no 
reason, a priori, to expect a large contribution to the optical 
emission from a 12.5 Gyr OSP. Therefore, in order to explore 
the possibility of YSPs dominating the optical light from 
the objects in our sample, we used a combination of two 
YSP components: a lYSP with ages in the range 0.3 - 2.0 
Gyr, and E{B - V) values of 0.0, 0.2 or 0.4, plus a VYSP 
with age in the range 1 - 100 Myr and reddenings (0.0 < 
E{B — V) < 2.0, increasing in steps of 0.1). Note that for 
this combination, we assume that the lYSPs have relatively 
low reddening {E(B - V) < 0.4). 
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The modelling results found for all the ULIRGs in the 
CS and the ES are shown in Tables [5l [71 [71 151 [9l and [TOl for 

Combination I, Combination II and Combination III (here- 
after Comb I, Comb II and Comb III). For the majority of 
the objects in our sample, the nebular correction mentioned 
before did not lead to a significant change in the modelling 
results for any of the combinations, and the uncertainties in 
the ages, reddenings and percentage contributions presented 
in the tables already account for such minor changes. There- 
fore we show only results for zero nebular correction in the 
tables for all objects with the exception of PKS1345-I-12. For 
this galaxy, the nebular correction changes the modelling re- 
sults when using Comb I for Ap A and to a lesser extent Ap 
B. In this case both the nebular corrected and uncorrected 
results are shown in Table [S] We now give a general overview 
of the results obtained for Comb I, Comb II and Comb III. 



IRAS 13539+2920 - 5kpc I: 12.5 Gyr OSP + 80Myr(E(B-V)=0.3) 




SoOO 3500 4000 4500 5000 5500 6000 6500 7000 

Wavelength (Angstroms) 



4 MODELLING RESULTS. 



4.1 Combination I 



In the case of Comb I, adequate fits (x^ < 1-0) to the over- 
all shape of the SED were obtained for all the extraction 
apertures modelled, with the exception of the 5 kpc aper- 
ture in IRAS 23060+0505, for which a power-law is required. 
However, although this combination fits the overall shape of 
the continuum, it fails to fit the absoption features (high- 
order Balmer Unes, Call K, G-band) for ~ 42% (54) of the 
extraction apertures modelled. Figure [3] shows a representa- 
tive example of the fits obtained using this combination. The 
aperture selected for the figure samples the south-eastern 
nucleus (5kpcl) of IRAS 13539-1-2920. The figure shows the 
best fitting model to both the continuum and the detailed 
absorption feartures, obtained for that aperture. It is clear 
from the figure that Comb I fits the optical continuum SED 
adequately, but there is an underprediction of the Call K 
absorption feature and, also, the fit to the Balmer lines and 
the G-band is not entirely satisfactory. This result suggests 
that the simple Comb I models may not be adequate for 
modelling the detailed absorption features for a subset of 
the ULIRGs in our sample. However, note that, the greatest 
disagreement between the Comb I fits and the data often oc- 
curs for the Call K line, which potentially has a contribution 
from absorption by the ISM in the host galaxjjfj. Therefore, 
for many of the 42% without adequate fits to the absorption 
lines, it is difficult to rule out the idea that Comb I provides 
adequate fits to the stellar continuum. 

Overall, even if the Comb I fits are not always perfect, 
they give single- valued estimates of the luminosity-weighted 
age, reddening and percentage contribution of the YSPs 
within each aperture. Therefore the results of Comb I - pre- 
sented in Table [5] and [7]- are still useful when performing a 
statistical analysis of the results. 



^ ~ 46% of the apertures for which Comb I does not perfectly fit 
the Call K absorption feature show an Nal D excess. The Nal D 
line is also potentially affected by the ISM of the host galaxy. 




700 3760 3600 3860 3900 3960 4000 4060 4100 4160 4200 4260 4300 

Wavelength (Angstroms) 




4200 4240 4280 4320 4360 4400 4440 4480 

Wavelength (Angstroms) 

Figure 3. An example of the fits obtained using Comb I for 
the south-western nucleus (aperture 5 kpc I) in the case of IRAS 
13539+2920. The model shown in the figure comprises a 12.5 Gyr 
OSP plus an 80 Myr VYSP {E(B - V) = 0.3) which contribute 
17% and 79% respectively to the flux in the normalising bin (4600 
- 4700 A). The green, blue and red spectra correspond to the OSP, 
the VYSP, and the sum of the two components respectively. It 
is clear from the figure that the shape of the continuum is ade- 
quately fitted. However the fit to the high order Balmer lines and 
the G-band is not entirely satisfactory (although some overpre- 
diction is expected for the Balmer lines, since the spectrum shown 
in the flgure has not been corrected for nebular emission) . More- 
over, there is a clear underprediction of the Call K absorption 
line. 



12 J.Rodriguez Zaurm, C.N.Tadhunter and R.M. Gonzalez Delgado^ 

IRAS 13539+2920 - 5kpc I: 500Myr(E(B-V)=0.2) + p.l. IRAS 13539+2920 - 5kpc I: 500Myr(E(B-V)=0.2)+6Myi(E(B-V)=0.5) 




3000 3500 4000 4500 6000 6500 6000 6600 7000 

Wavelength (Angstroms) 




3000 3500 4000 4500 5000 6500 6000 6500 7000 

Wavelength (Angstroms) 
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4200 4240 4280 4320 4360 4400 4440 4480 
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Figure 4. Same as Figure |3] using Comb II. The figure shows a 
model comprising a 500 Myr lYSP {E(B — V) = 0.2) and a power- 
law (power-law index a = -1.06) which contribute 78%, 21% re- 
spectively to the fiux in the normalising bin (4600 - 4700 A), with 
no need of a contribution from an 12.5 Gyr OSP. The green, blue 
and red spectra correspond to the power-law, the lYSP, and the 
sum of the two components respectively. It is clear from the figure 
that the underprediction of the Call K absorption line disppears 
when a power-law is included. A small improvement in the fit to 
the G-band is also present using this combination. 



a 3 




4200 4240 4280 4320 4360 4400 4440 4480 

Wavelength (Angstroms) 

Figure 5. Same as Figure [3] but using Comb III. The model 
shown in the figure is comprised of 500 Myr lYSP {E{B — V) = 
0.2) plus a 6 Myr VYSP {E{B-V) = 0.5), which contribute 56% 
and 40% respectively to the flux in the normalising bin (4600 - 
4700 A). The green, blue and red spectra correspond to the VYSP, 
the lYSP, and the sum of the two components respectively. In the 
case of Comb III, both the continuum and the absorption features 
(high-order Balmer lines. Call K and the G-band) are better fitted 
using this combination than in the case of Comb I and Comb II. 
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4.2 Combination II 

For those ULIRGs in our sample also classified as Sy2 galax- 
ies, it is possible that the optical spectrum is a compos- 
ite of stellar and AGN continuum emission (e.g. scattere d 
light or weak direct AGN light: see Tadhunter et al.l , [2002h . 
In this context, we note that significant starburst activity 
has alrea dy been detected in th e nucl ear regions in Sy2 
galaxies. iGonzalez Delgado et al.l |20o3) found high-order 
Balmer lines and He I absorption lines, strong indicators of 
the presence of YSPs, in their optical spectroscopic study 
of a sample of 20 Sy2 galaxies. Most of the objects in 
the lGonzalez Delgado et al.l (|200ll ) sample are classified as 
LIRGs, and one, Mrk 273, is classified as a ULIRG and in- 
cluded in the work presented here. Therefore, it is straight- 
forward to understand that Comb II, including a power law, 
is particularly useful for those ULIRGs in our sample also 
classified as Sy2 galaxies at optical wavelegths. 

Since the power law also represents a highly reddened 
VYSP, comb II gives an idea of the contribution level of the 
OSPs in each galaxy. Therefore, we decided to use it for the 
whole sample of ULIRGs. Adequate fits (x^ < 1.0) to the 
overall SED shapes were found for all of the apertures of all 
objects using this combination. On the other hand, this com- 
bination fails to fit the detailed absorption features for 10% 
(13 apertures) of the apertures modelled. Figure |3] shows an 
example of a fit obtained for IRAS 13539-1-2920 using Comb 
II for the same extraction aperture shown in Figure |3l It is 
clear from the figure that the underprediction of the Call K 
absorption line disappears when a power-law is included. Fi- 
nally, as shown in Tables [7] and [8] the minimum precentage 
contribution of the OSP component is less that 10% (due to 
the uncertainties inherent in the modelling techniques per- 
centage contributions < 10% are considered negligible) for 
most of the apertures. This result suggest that an OSP com- 
ponent is not required to fit the data in the majority of the 
objects in the ES. 



4.3 Combination III 

A drawback of Comb II is that, for those cases in which 
the power-law is likely to represent a VYSP component, it 
gives no information about the detailed properties (age and 
reddening) of such a component. Therefore, in order to bet- 
ter investigate the detailed properties of the VYSPs in the 
galaxies, we used Comb III. Adequate fits to both the con- 
tinuum and the detailed absorption features are obtained 
for all but 9 (~ 7%) of the extraction apertures modelled. 
In the minority of the cases for which no adequate fits could 
be obtained using Comb III, either an OSP or a power- 
law component was required in addition to a YSP. Figure [S] 
shows an example of a fit obtained using Comb III for the 
same extraction aperture shown in Figures [3] and |4] in the 
case of IRAS 13539+2920. The figure clearly demonstrates 
that Comb III provides the best fit to all of the detailed 
absorption features. 

At this stage, it is worth re-iterating that there is no 
reason, a priori, to expect a major contribution to the optical 
emission from a 12.5 Gyr OSP in ULIRGs. The lack of an 
OSP is also suggested by the results of Combination II and, 
to a lesser extent, those of Combination I. Therefore, in the 
following, we will concentrate on the modelling results for 



Comb III and use those for Comb I and II as additional 
information, unless otherwise specified. 



5 ESTIMATING MASSES AND BOLOMETRIC 
LUMINOSITIES 

It is also possible to estimate the masses and bolometric lu- 
minosities associated with the stellar populations detected 
at optical wavelengths using the modelling results. For this 
purpose, we extracted and modelled a large aperture includ- 
ing all the structure visible in the 2-D frames (and in the 
spatial cut in Figure [2]) for each object. Assuming that the 
stellar population mix remains constant across the entire 
visible extents of the galaxies, we chose a model from all 
the adequate fits obtained, selected to be as representative 
as possible of all the stellar populations in the galaxy. Once 
the model was selected, we estimated the mass contribu- 
tion of each of the stellar components in the large aperture. 
Since the slit doe s not cover the fu ll extent of each galaxy, 
we then used the lKim et al.l (2002) total magnitudes of the 
galaxies in the R-band to scale the results and estimate the 
masses in the various stellar components. However, due to 
the diversity of the stellar populations found, it is some- 
times hard to select a particular representative model for 
some ULIRGs. In such cases, the selection is biased towards 
models consistent with the results obtained for the nuclear 
region, since it is likely that most of the mid- to far-IR lumi- 
nostity arises from these regions. The estimated total masses 
of the YSPs (VYSPs plus lYSPs, resuhs presented in Table 
10) faU in the range 0.18x10^° < Mysp < 50 xlO^° Mq. 
In the cases where both VYSPs and lYSPs are present, the 
lYSP is, in most cases, more massive than the VYSP (the 
exceptions are IRAS 10494-^4424, IRAS 14252-15550 and 
IRAS 17028+5817W) 

To estimate the bolometric luminosities associated with 
the stellar populations detected in the optical, we assumed 
the same combination of stellar populations as used for 
the mass measurements, scaling the unreddened stellar tem- 
plates using the mass values previously obtained. We then 
integrated the scaled synthetic spectra over their entire 
wavelength range. The results are presented in Tabled We 
find stellar bolometric luminosities in the range of 0.07 x 
10^^ < Lboi < 2.2 X 10^^ Lq. Assuming that most of the 
optical light is absorbed and then re-processed by dust, we 
can compare the estimated bolometric luminosities with the 
mid- to far-IR luminosities of the sources, presented in Ta- 
ble [1] Column 7 in the table shows the estimated bolometric 
luminosities of the YSPs detected in the optical expressed 
as percentages of the measured mid- to far-IR luminosities. 

Finally, we would like to emphasize that the results in 
Table |3] are indicative rather than exact, and they must 
be interpreted carefully due to the various uncertainties at- 
tached to the process. Above all, the major source of un- 
certainity arises from the selection of the model; different 
combinations of stellar populations may change the results 
significantly. 
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54 
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Table 4. Estimated total masses and bolometric luminosities associated with the stellar components detected at optical wavelengths. 
Col (1): object name. Col (2): estimated mass associated with the VYSP component. Col (3): estimated masses associated with the 
lYSP component. Col (4): estimated masses associated with the OSP component (if present). Col (5): es timated bolometric luminosities 
associated with the stellar components. Col(6): mid- to far IR luminosity of the sources from lKim fc San ders (199^. Col (7): estimated 
bolometric luminosities of the YSPs detected in the optical expressed as percentages of the mid- to far-IR luminosities. 
Note that since it was not possible to model the stellar populations in the nuclear regions of the three ULIRGs classified as Syl galaxies 
(IRAS 12540+5708,IRA 15462-0450 and IRAS 21219-1757), no results for these objects are presented in the table (see Apendix A for 
details). 



6 SUMMARY OF GENERAL PATTERNS AND 
RESULTS 

In order to give an overview of the properties of the stellar 
populations of the ULIRGs in our sample as a whole, the 
general results and patterns obtained from the modelling 
can be summarized as follows. 

• Modelling technique: three different combinations of 
stellar populations have been used for the analysis presented 
here: combination I (OSP + YSP), combination II (OSP 
+ YSP + p.l.) and combination III (VYSP + lYSP). In 
general, adequate fits to the overall shapes of the continua 
in individual apertures are obtained using all of the combi- 
nations. However, Comb I fails to fit the detailed absorption 
features for 42% (54) of the extraction apertures. This per- 



centage is reduced down to 10% (13 apertures) and 7% (9 
apertures) in the cases of Comb II and Comb III, respec- 
tively. Note that adequate fits for both the continuum and 
the detailed absorption features are found for 60% of the ex- 
traction apertures using alloi the combinations. This clearly 
demonstrates the difficulty in obtaining a unique solution 
when using combinations of different stellar populations to 
model the optical spectra of starburst systems. 

• Young stellar populations: young stellar popula- 
tions (YSPs, tvsp < 2.0 Gyr) are present in all apertures in 
all sources, with the exception of IRAS 23327+2913 5kpc- 
II. Furthermore, very young stellar populations (VYSPs, 
tvYSP < 100 Myr) are essential to obtain adequate fits in 
all hut 20 (~ 15%) of the extraction apertures, 14 of which 
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sample the extended regions of Arp220 and the radio galaxy 
IRAS 13451+1232 (PKS1345+12). 

At this stage we would like to add a caveat about the prop- 
erties of the so-called intermediate-age young stellar popu- 
lations (lYSPs, 0.1 < tiYSP < 2.0 Gyr). For this component, 
we have assumed a limited range of ages and, in particular, 
reddening values, in comparison with those of the VYSPs. 
Therefore, the ranges of age and reddening presented in Ta- 
bles [5] and [TU] for the lYSPs are partially a consequence of 
these assumptions. 

• Old stellar populations: old stellar populations 
(OSPs, 12.5 Gyr) do not make a major contribution to the 
optical light in most apertures for most sources. The mod- 
elling results show that OSPs are essential for fitting the 
optical continuum in only 7 ('^ 5%) of all the apertures con- 
sidered; in the majority of the cases, adequate fits can be 
obtained using combinations of lYSPs and VYSPs. Only in 
the cases of IRAS 23327-f 2913 (5kpc-II) and IRAS 21208- 
0519 (5kpc-II) we find that the OSPs dominate. In addition, 
in the cases of IRAS 13451-^1232, IRAS 14394-^5332 5kpc- 
II, IRAS 16487-h5447E and IRAS 23234-^0946 5kpc-II either 
OSPs, or dominant "old" lYSPs (1-2 Gyr), are required to 
model the data. 

• Reddening: the study presented here further empha- 
sizes the importance of carefully accounting for dust redden- 
ing effects when modelling the stellar populations of star- 
forming galaxies. 

• Double nucleus systems: 20 of the 36 ULIRGs in 
our sample are classified as double nucleus systems and one, 
IRAS 14394+5322, is classified as a multiple (> 2 nucleus) 
system. For 14 of these systems, it was possible to study the 
stellar populations within the different nuclei individually. 
In 71% of the cases (10 of the 14 objects) the results are 
consistent with the stellar populations being similar in the 
two nuclei in terms of ages, reddenings and percentage con- 
tributions. However, in four cases (IRAS 14394+5332, IRAS 
16487+5447, IRAS 21208-0519 and IRAS 23327+2913) we 
find marked differences between the mix of ages of the stellar 
populations in the two nuclei. 

• Extended and nuclear apertures: in the majority 
of the objects we find that we require the same stellar com- 
ponents (lYSP and VYSP) at all locations of the systems to 
model the spectra adequatelly. 

• Masses of the YSP: the estimated total masses of the 
YSPs (VYSPs plus lYSPs) faU in the range 0.18x10^° < 
Mysp < 50 xlO^° M0. In the cases where both VYSPs and 
lYSPs are present, the lYSP is, in most cases, more massive 
than the VYSP (the exceptions are IRAS 10494+4424, IRAS 
14252-15550 and IRAS 17028+5817W). 

• Bolometric luminosities: we have also estimated the 
bolometric luminosities associated with these stellar popu- 
lations. For 48% of the objects in the ES sample (16 of the 
33 objects, excluding the Syl galaxies), the bolometric lu- 
minosities of the YSPs detected in the optical can account 
for a large fraction ( > 50%) of the mid- to far-IR luminosity 
of the source. 

More detailed information about individual sources is pre- 
sented in Appendix A. 



7 CONCLUSIONS 

This is the first of two papers presenting a detailed study of 
the stellar populations in a large sample of ULIRGs, based 
on long-slit optical spectroscopic observations. In this paper 
we have introduced the sample and described the analysis 
techniques and the general results obtained from the mod- 
elling. In addition, we present estimated values for the total 
masses of the stellar populations within the objects and the 
bolometric luminosities associated with them. 

Overall, the results for the sample in general are consis- 
tent with those obtained in the detailed studies of RZ07 and 
RZ08 for PKS1345+12 and Arp 220 respectively, in terms 
of: 

• the requirement for a mixture of VYSPs and lYSPs, 
with the proportional contribution of the VYSP component 
increasing towards the nuclei in most cases; and 

• the importance of taking into account the reddening 
of the YSP component, with clear evidence of increasing 
reddening towards the nuclei in a number of the objects. 

The work presented here clearly shows the utility of 
spectral synthesis modelling techniques for studying the 
properties of the stellar populations in star forming galax- 
ies. In addition, it also shows the difficulty in obtaining a 
unique solution when using combinations of different stellar 
populations to model spectra, especially on the basis of fits 
to the SED alone (i.e. not examining the detailed absorption 
features). We have shown that some results, and therefore 
the conclusions based on them, can change significantly if 
different models are used. 

The results described in this paper will be analyzed in 
more detail in a forthcoming paper (Rodriguez Zaurm et 
al. 2009b), where they will be discussed in the context of 
evolutionary models for merging systems, and the results 
obtained for similar systems in the high redshift universe. 
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Figure 6. The extracted ID-spectra for all the ULIRGs and apertures in the ES. When the dichroic effects are important, it is 
not possible to show an adequate plot of the entire wavelength range, and the dichroic region is removed from the plots. The fluxes 
are presented in wavelength units. 
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Table 5. Modelling results for the CS using Comb I (12.5 Gyr OSP + YSP). Col (1): object name. Col (2): aperture label as indicated in 
Figure |2l Col (3): this column indicates whether or not Comb I fits both the overall shape of the continuum and the detailed absorption 
features: Y = yes, N = No. Col (4): range of ages for the YSP component that adequately fit the data (x^ < 1.0). Col (5): range of 
E(B — V) values of the YSPs in Col (4). Col (6): the percentage flux contribution of the YSP component to the model in the normalising 
bin. 
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E(B-V) 


%YSP 


Name 




features 


YSP 










IRAS 




fitted? 


(Gyr) 










(1) 


(2) 


(3) 


(4) 




(5) 




(6) 


08572+3915 


5kpc 


N 


0.06 - 0.1 


0.0 - 


0.3 


65 


- 90 




Ap A 


Y 


0.1 - 0.3 


0.0 - 


0.1 


45 


- 55 




Ap B 


N 


0.07 - 0.1 


0.0 - 


0.2 


60 


- 85 




Ap C 


Y 


0.07 - 0.1 


0.2 - 


0.3 


90 - 


100 


10190+1322 


5kpc-I 


N 


0.05 - 0.08 


0.3 - 


0.5 


65 


- 90 




5kpc-II 


N 


0.2 




0.5 




85 




Ap A 


Y 


0.4 - 0.7 


0.0 - 


0.2 


60 


- 95 




Ap B 


N 


0.04 - 0.06 


0.5 - 


0.6 


80 


- 95 


10494+4424 


5kpc 


N 


0.07 - 0.1 




0.7 


95 - 


100 




Ap A 


Y 


0.4 


0.0 - 


0.1 


72 


- 87 




Ap B 


N 


0.04 - 0.1 


0.0 - 


0.4 


75 - 


100 


13305-1739 


5kpc 


N 


0.05 - 0.1 


0.0 - 


0.4 


65 


- 90 


14252-1550 


5kpc 


N 


0.05 - 0.1 


0.5 - 


0.6 


80 


- 85 




Ap A 


Y 


0.2 - 0.7 


0.0 - 


0.2 


50 - 


100 


16156+0146 


5kpc-I 


Y 


0.01 




0.5 




75 




5kpc-II 


N 


0.08 - 0.1 


0.1 - 


0.3 


75 


- 98 




Ap A 


Y 


0.004 - 0.007 


0.3 - 


0.5 


40 


- 80 








0.03 - 0.04 


0.1 - 


0.3 


50 


- 75 


17044+6720 


5kpc 


Y 


0.05 - 0.08 




0.4 


80 


- 85 




Ap A 


Y 


0.05 - 0.1 


0.0 - 


0.3 


70 


- 90 




Ap B 


N 


0.05 - 0.1 


0.0 - 


0.2 


77 


- 95 


17179+5444 


5kpc 


N 


0.07 - 0.1 


0.4 - 


0.6 


50 


- 85 


23060+0505 


5kpc 


N 














Ap A 


Y 


0.05 - 0.1 


0.1 - 


0.4 


43 


- 65 




Ap B 


Y 


0.04 - 0.06 


0.2 - 


0.3 




65 


23389+0303 


5kpc I 


N 


0.004 




0.9 




-100 








0.04 




0.6 




90 




5kpc II 


N 


0.07 - 0.1 


0.2 - 


0.4 


70 


- 85 




Ap A 


N 


0.06 




0.5 




85 




Ap B 


N 


0.005 - 0.006 


0.8 - 


0.9 




85 




Ap C 


N 


0.05 - 0.08 




0.4 


63 


-76 




Ap D 


N 


0.1 




0.2 




72 



Table 6. Same as Table[5]but for the 10 additional ULIRGs included in the ES. 
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Object Name 


Aperture 


Age 


of YSP 


E(B 


- V) 


%YSP 


%p.l 


%OSP 


IRAS 




(Gyr) 


















(1) 


(2) 


(3) 




(4) 




(5) 




(6) 




(7) 




5kpc 


0.3 - 


0.6 


0.0 - 


0.2 


48 - 


60 


20 - 


35 


5 - 


30 




Ap A 


0.4 - 


1.0 


0.0 - 


0.3 


50 - : 


LOO 


6 - 


12 


< 


25 




Ap B 


0.2 - 


0.6 


0.0 - 


0.2 


66 - 


75 


10 - 


25 


< 


24 


nni 88-0856 




0.3 - 


0.6 


0.3 - 


0.4 


75 - 


90 


8 - 


10 


< 


16 




5kpc 


0.03 


- 0.1 


0.0 - 


0.5 


30 - 


50 




-40 


10 - 


30 


12072-0444 


Skpc"'^ 


- 










- 




- 




- 




Ap A 


0.2 - 


0.6 


0.0 - 


0.2 


40 - 


65 


7 - 


30 


< 


50 




Ap B 


0.5 - 


0.6 


0. 





50 - 


60 




40 


< 


10 


1211 2-i-n'^n5 

-L ^ -L -L \ W O W ^J 




0.4 - 


0.5 


0.0 - 


0.1 


50 - 


60 




-40 


< 


10 




5kpc II 


0.2 - 


0.5 


0.0 - 


0.2 


50 - 


60 


30 - 


40 


< 


15 




Ap A 


0.2 - 


0.8 


0.0 - 


0.1 


40 - 


70 


20 - 


40 


< 


45 


1 2540-1-5708 




- 










- 




- 




- 




Ap A 


0.4 - 


1.4 


0.1 - 


0.5 


40 


-75 


25 - 


40 


< 


35 




Ap B 


0.2 




0. 







40 




30 




30 




Ap C 


0.09 


- 0.2 


0.0 - 


0.1 


70 


-88 


10 - 


30 









Ap D 


- 










" 




" 




" 




Ap E 


0.1 - 


0.2 


0.0 - 


0.1 




30 




40 




30 


1 '^428-1-5608 




0.5 - 


0.8 


0.0 - 


0.1 


45 - 


60 




-40 


< 


15 




Ap A 


0.4 - 


0.8 


0.0 - 


0.3 


50 - 


85 


10 - 


15 


< 


35 




Ap B 


0.2 - 


0.7 


0.0 - 


0.2 


33 - 


64 


6 - 


15 


30 - 


65 




Ap C 


0.4 - 


0.7 


0.0 - 


0.2 


50 - 


90 


6 


- 9 


< 


45 




Ap D 


0.5 - 


0.8 


0.0 - 


0.3 


50 - 


80 


18 - 


23 


< 


30 




Ap E 


0.5 - 


0.8 


0.0 - 


0.1 


40- 


60 




-40 


< 


20 




Ap F 


0.5 - 


0.7 


0.0 - 


0.4 


45 - 


80 


15 - 


30 


< 


30 


1345H-1232(PA160) 


5kpc'^ 












- 




- 




- 




Ap A 


0.2 - 


0.8 


0.0 - 


0.2 


20 - 


35 




- 5 


60 - 


75 




Ap B^ 












- 




- 




- 




Ap C 


0.4 - 


0.8 


0.0 - 


0.3 


35 - 


75 


15- 


25 


5 - 


45 




Ap D 


0.4 - 


1.2 


0.0 - 


0.5 


30 - 


85 


10 - 


16 


< 


45 


1 3451 -1-1 232('PA230"l 


Ap A 


0.5 - 


1.4 


0.0 - 


0.5 


40 - 


90 


3 - 


10 


< 


50 




Ap B 


0.4 - 


1.0 


0.0 - 


0.2 


25 - 


50 


5 


-10 


45 - 


65 


1 353Q4-2Q20 




0.5 - 


0.6 


0.0 - 


0.2 


65 - 


80 


20 - 


30 









5kpc-II 


0.5 - 


0.8 


0.2 - 


0.5 


50 - 


70 


35 - 


40 


< 


15 


14060-1-2919 




0.004 


0. 


5 




93 




6 









Ap A 


0.3 - 


0.8 


0.0 - 


0.4 


60 - 


90 


10 - 


15 


< 


25 




Ap B 


0.04- 


0.06 


0.0 - 


0.4 


60 - 


85 


10 - 


25 


< 


10 




Ap C 


0.2 - 


0.3 


0. 







50 




40 




11 


14S48-1447 


5kpc-I^ 
























5kpc-II 


0.5 




0. 







62 




38 









Ap A 


0.1 - 


0.3 


0.0 - 


0.1 


40 - 


50 




40 


10 - 


20 




Ap B 


0.2 - 


0.6 


0.0 - 


0.1 


30 - 


40 


30 - 


40 


12 - 


40 




Ap C 


0.09 


- 0.4 


0.0 - 


0.4 


45 - 


80 


5- 


35 


< 


35 




Ap D 


0.05 


- 0.4 


0.0 - 


0.5 


40 - 


75 


5 - 


30 


< 


50 


14394-1-5332 


5kpc-li 
























5kpc-II 


0.7 - 


2.0 


0.0 - 


0.2 


40 - 


70 


30 - 


35 


< 


30 




Ap A 


0.3 - 


0.8 


0.0 - 


0.2 


42 - 


60 


30 - 


35 


< 


25 




Ap B 


0.4 - 


1.0 


0.0 - 


0.3 


40 - 


60 


23 - 


35 


< 


30 


15130-1-1958 


5kpc 


0.2 - 


0.5 


0.0 - 


0.2 


55 - 


70 




30 


< 


15 


15206-f3342 


5kpc'^ 























Table 7. Modelling results for the CS using Comb II (OSP + YSP -I- power-law). Col (1): object name. Col (2): aperture label as 
indicated in Figure [2] Col (3): range of ages for the YSP components that adequately fit the data (both the continuum and the detailed 
features). Col (4): range of E{B — V) values of the YSPs in Col (3). Col (5): the percentage flux contribution of the YSP component to 
the model, in the normalising bin. Cols (6) and (7): same as Col (5) for the power-law and OSP respectively. 
^ No adequate fits (fitting both the continuum and the detailed absorption features) are found using Comb II. 

^ No adequate fits are found with a non negligible contribution of the p.l. component (i.e., different than those of Comb I) and therefore, 
no results are presented in the table. 

^ When a p.l. is included, adequate fits are found for the entire range of YSP combinations used and, therefore, no results are presented 
in the table. 
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Object Name 


Aperture 


Age 


of YSP 


E(B - V) 


%YSP 


%p.l 


%OSP 


IRAS 




(Gyr) 














(1) 


(2) 


(3) 




(4) 




(5) 




(6) 


(7) 


153274-234nCPA160'l 


5kpc 


0.4 - 


0.7 


0.4 - 0.5 


64 - 


85 


15 - 


18 


<20 




Ap A 


0.5 - 


0.7 


0.0 - 0.1 


75 - 


98 




~2 


<20 




Ap B 


0.5 - 


0.7 


0.1 - 0.2 


85 - 


95 


3 


- 5 


<20 




Ap C 


0.4 - 


0.7 


0.3 - 0.5 


70 - 


98 


4 


- 5 


<20 




Ap D 


0.4 - 


0.6 


0.5 - 0.8 


62 - 


75 


27 - 


35 


<8 




Ap E 


0.4 - 


0.6 


0.6 - 0.8 


75 - 


85 


10 - 


16 


<15 




Ap F 


0.1 - 


0.7 


0.2 - 0.4 


53 - 


85 


4 - 


17 


<43 




Ap G 


0.4 - 


0.7 


0.0 - 0.2 


55 - 


90 


4 


- 8 


<40 




ApxOTAL 


0.3 - 


0.7 


0.2 - 0.4 


48 - 


87 


5 - 


15 


<45 




Ap A 


0.3 - 


0.7 


0.1 - 0.3 


60 - 


93 


3 


- 7 


<37 




Ap B 


0.2 - 


0.6 


0.2 - 0.5 


50 - 


75 


6 - 


18 


10 - 44 




Ap C 


0.3 - 


0.6 


0.4 - 0.5 


65 - 


85 


6 - 


12 


8-27 




Ap D 


0.3 - 


0.7 


0.2 - 0.4 


61 - 


95 


3 


- 9 


<37 




Ap E 


0.3 - 


0.7 


0.1 - 0.4 


52 - 


85 


4 - 


10 


9-45 




Ap F 


0.3 - 


0.9 


0.1 - 0.5 


30 - 


85 


2 - 


13 


10 - 70 




Ap G 


0.2 - 


1.0 


0.0 - 0.3 


38 - 


90 


7 - 


16 


<60 




A nT^/~\T^ A T 

1 1 ALi 


0.2 - 


0.6 


0.2 - 0.4 


57 - 


75 


4 - 


12 


12 - 40 


153274-234f)CPA75*'l 


Ap A 


0.3 - 


0.6 


0.0 - 0.2 


55 - 


90 


3 - 


10 


<40 




Ap B 


0.4 - 


0.7 


0.0 - 0.2 


59 - 


95 


2 


- 5 


<41 




Ap C 


0.4 - 


0.6 


0.0 - 0.2 


63 - 


90 


4 


- 8 


<33 




Ap D 


0.4 - 


0.7 


0.0 - 0.2 


55 - 


90 


7 - 


10 


<40 




Ap E 


0.4 - 


1.2 


0.0 - 0.3 


40 - 


90 


9 - 


11 


<52 




Ap F 


0.5 - 


1.6 


0.0 - 0.3 


45 - 


95 


2 - 


10 


<48 




ApxOTAL 


0.4 - 


0.7 


0.0 - 0.3 


51 - 


93 


7 


- 9 


<45 




5kpc-I^ 


- 








- 




- 








0.2 - 


0.4 


00 - 0.2 


50 - 


65 


20 - 


30 


< 25 




Ap A^ 


- 








- 




- 




Ifi487+5447NE 




0.5 - 


1.0 


0.1 - 0.3 


30 - 


50 


10 - 


15 


40 - 60 




Ap A 


0.04 


- 0.1 


0.0 - 0.3 


25 - 


45 


15 - 


33 


33 - 40 


1 fi4S74-'i447NW 


5kpc"^ 


0.04 




0.2 




60 




30 


10 




Ap A 


0.2 - 


0.4 


0.0 - 0.1 


48 - 


55 


16 - 


20 


25 - 30 




Ap B 


0.04 


- 0.3 


0.0 - 0.2 


50 - 


66 


10 - 


35 


20 - 40 


1 7028+581 7NR 


5kpc^ 


- 








- 




- 




17028+5817NW 


5kpc 


0.3 - 


0.6 


0.1 - 0.4 


50 - 


60 


25 - 


35 


10 - 25 




Ap A 


0.3 - 


0.6 


0.0 


64 - 


75 


16 - 


20 


< 15 


20414-1651 


5kpc 


0.3 - 


0.6 


0.0 - 0.2 


50 - 


65 


20 - 


40 


< 15 


21 208-051 Q 




0.1 - 


0.2 


0.0 - 0.2 


45 - 


50 


30 - 


40 


10 - 25 




Sknr-II 


0.2 - 


1.2 


0.0 - 0.5 


25 - 


85 


10 - 


22 


< 65 




Ap A 


0.2 - 


0.6 


0.0 - 0.2 


46 - 


70 


5 - 


20 


10 - 50 




Ap B 


0.2 - 


0.6 


0.0 - 0.1 


47 - 


55 


25 - 


35 


10 - 22 


22491-1808 


5kpc^ 










- 




- 






An A 


0.4 - 


0.6 


no 9 


67 - 


85 


15 - 


30 


n 
u 




Ap B^ 










- 




- 






Ap 




















5kpc 


0.05 


- 0.5 


no n 
u.u - u.o 


25 - 


45 


15 - 


40 






Ap A 




















Ap B 


0.07 


- 0.4 


0.0 - 0.3 


35 - 


44 


6 - 


20 


40 - 50 


23234+0946 


5kpc-I 


0.04 


- 0.2 


0.0 - 0.4 


20 - 


63 


17- 


40 


20 - 40 




5kpc-II 


0.1 - 


1.0 


0.0 - 0.4 


35 - 


80 


5 - 


20 


< 60 




Ap A 


0.5 - 


0.8 


0.0 


65 - 


85 


10 - 


20 


< 15 


23327+2913 


5kpc-I 


0.001 - 0.003 


0.5 


40 - 


45 


25 - 


35 


25 






0.01 




0.1 




30 




40 


30 




5kpc-Il3 




















Ap A 


0.001 - 0.003 


0.5 


40 - 


45 


25 - 


35 


25 






0.01 




0.1 




30 




40 


30 




Ap B3 



















Table 7. Continued 
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Object Name 


Aperture 


Age of YSP 


E(B - V) 


%YSP 


%p.l 


%OSP 


IRAS 




(Gyr) 














(1) 


(2) 


(3) 


(4) 




(5) 




(6) 


(7) 


08572+3915 


5kpc 


0.2 - 0.6 


0.0 - 0.2 


60 


- 70 


30 


- 36 


< 25 




Ap A 


0.2 - 0.6 


0.0 - 0.3 


50 


- 85 


5 


- 20 


< 45 




Ap B 


0.4 


0.0 - 0.1 


65 


- 75 


25 


- 30 







Ap C 


0.09 - 0.4 


0.0 - 0.3 


73 


- 95 


6 


- 35 


< 12 


1 01 QO-l-1 322 

-L U -L W 1 _1_ ^ ^ 


5kpc-I^ 
5kpc-II^ 








- 
- 




- 
- 






Ap A 


0.4 - 0.7 


0.0 - 0.2 


70 


- 95 




~ 5 


< 25 




Ap 








- 




- 




1 n4Q4-l-4424 


Ap A^ 








- 

- 




- 

- 






Ap B 


0.07 - 0.2 


0.0 - 0.2 


65 


- 85 


10 


- 40 


< 10 


13305-1739 


5kpc 


0.3 - 0.5 


0.0 - 0.2 


50 


- 60 




40 


< 10 


14252-1550 




0.2 - 0.6 


0.0 - 0.2 


40 


- 66 


15 


- 35 


< 35 




Ap A 


0.2 - 0.8 


0.0 - 0.2 


60 


- 85 


10 


- 20 


< 45 


16156+0146 




< 0.03 


.2 - 0.7 


25 


- 55 


30 


- 40 


14 - 40 




5kpc-II 


0.4 - 0.5 


0.0 - 0.1 


70 


- 80 


20 


- 26 







Ap A 


0.05 - 1.0 


0.2 - 0.5 


15 


- 66 


30 


- 40 


< 30 


17044+6720 


5kpc 


0.2 - 0.4 


0.0 - 0.1 




- 50 


36 


- 40 


10 - 15 




Ap A 


0.2 - 0.7 


0.0 - 0.2 


50 


- 62 


35 


- 40 


< 15 




Ap B 


0.4 - 0.5 


0.0 




55 




40 


5 


17179+5444 


5kpc 


0.2 - 0.6 


0.1 - 0.3 


40 


- 62 


22 


- 35 


6 - 50 


23060-1-0505 


5kpc 


0.01 


0.4 




72 




6 


22 




Ap A 


0.07 - 0.6 


0.0 - 0.3 


35 


- 50 


6 


- 30 


20 - 35 




Ap B 


0.05 - 0.6 


0.0 - 0.4 


40 


- 60 


8 


- 35 


5 - 30 


23389+0303 


5kpc-li 


















5kpc-II 


0.2 - 0.7 


0.0 


53 


- 70 


12 


- 30 


< 35 




Ap A 


0.4 - 0.6 


0.0 - 0.1 




- 60 




- 40 







Ap Bi 


















Ap C 


0.5 - 0.8 


0.0 - 0.1 


50 


- 70 




- 30 


<20 




Ap D 


0.2 - 0.5 


0.0 - 0.1 


55 


- 65 


10 


- 30 


2 - 30 



Table 8. Same as Table [7] but for the 10 additional ULIRGs included in the ES. 
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Object 


Aperture 


Age of 


E{B - V) 


Age of 


E(B 


-V) 


%VYSP 


Name 




lYSP 


of 


VYSP 




of 






IRAS 




(Gyr) 


lYSP 


(Myr) 


VYSP 






(1) 


(2) 


(3) 




(4) 


(5) 




(6) 




(7) 


00091-0738 


5kpc 


0.5 - 


1.0 


0.0 - 0.4 


6 - 50 


0.5 - 


1.0 


30 - 


65 




Ap A 


0.7 - 


2.0 


0.0 - 0.2 


< 100 


0.0 - 


1.5 


< 


50 




Ap B 


0.3 - 


2.0 


0.0 - 0.4 


< 100 


0.2 - 


1.0 


10 - 


90 


00188-0856 


5kpc 


0.5 - 


0.7 


0.0 - 0.4 


< 100 


0.2 - 


1.6 


< 


40 


01004-2237 


5kpc 


0.5 - 


0.7 


0.2 - 0.4 


9 - 20 


0.1 - 


0.3 


60 - 


75 


12072-0444 


5kpc 


0.7 - 


2.0 


0.0 - 0.4 


7 - 20 


0.4 - 


0.6 


45 - 


65 




Ap A 


0.5 - 


2.0 


0.0 - 0.4 


< 100 


0.0 - 


1.2 


20 - 


65 




Ap B 


1.0 - 


2.0 


0.0 


40 - 100 


0.0 - 


0.3 


40 - 


60 


12112-1-0305 


5kpc-I 


0.5 - 


1.0 


0.0 - 0.4 


7 - 40 


0.5 - 


0.8 


50 - 


90 




5kpc-II 


0.3 - 


2.0 


0.0 - 0.4 


4 - 60 


0.0 - 


0.5 


35 - 


90 




Ap A 


0.5 - 


2.0 


0.0 - 0.4 


< 50 


0.6 - 


1.0 


35 - 


80 


12540-1-5708 


5kpc 


- 




- 


- 




- 




- 




Ap A 


2 




0.0 


5 - 10 


0.0 - 


0.1 


30 - 


35 




Ap B 


0.5 - 


2.0 


0.0 - 0.4 


10 - 100 


0.1 - 


0.5 


30 - 


65 




Ap C 


0.3 - 


2.0 


0.0 - 0.4 


50 - 80 


0.0 - 


0.1 


70 - 


96 




Ap D 


- 




- 


- 




- 




- 




Ap E 


0.7 - 


2.0 


0.0 - 0.4 


6 - 70 


0.5 - 


1.0 


40 - 


75 


13428-1-5608 


5kpc 


0.7 - 


2.0 


0.0 - 0.4 


6 - 50 


0.4 - 


0.8 


20 - 


50 




Ap A 


0.7 - 


2.0 


0.0 - 0.2 


< 100 


0.0 - 


0.6 


5 - 


50 




Ap B 


2.0 




0.0 


50 - 100 


0.4 - 


0.5 


25 - 


33 




Ap C 


0.7 - 


2.0 


0.0 - 0.2 


< 100 


0.0 - 


0.8 


10 - 


50 




Ap D 


0.7 - 


2.0 


0.0 - 0.4 


5 - 60 


0.3 - 


0.8 


25 - 


60 




Ap E 


0.7 - 


2.0 


0.0 - 0.4 


6 - 50 


0.4 - 


0.8 


20 - 


50 




Ap F 


0.7 - 


2.0 


0.0 - 0.4 


5 - 60 


0.3 - 


0.6 


25 - 


55 


13451-I-1232(PA160) 


5kpc 


2.0 




0.0 - 0.2 


< 100 


0.0 - 


2.0 


< 


10 




Ap A 


2.0 




0.0 - 0.1 


< 100 


0.0 - 


2.0 


< 


10 




Ap B 


2.0 




0.1 - 0.2 


< 100 


0.0 - 


1.4 


5 - 


20 




Ap C 


1.0 - 


2.0 


0.1 - 0.2 


< 100 


0.1 - 


0.8 


20 - 


35 




Ap D 


1.0 - 


2.0 


0.0 - 0.2 


< 100 


0.0 - 


0.9 


5 - 


15 


13451-I-1232(PA230) 


Ap A 


2.0 




0.0 - 0.1 


< 100 


0.0 - 


2.0 


< 


12 




Ap B 


2.0 




0.1 


< 100 


0.0 - 


2.0 


< 


11 


13539-f2920 


5kpc-I 


0.5 - 


2.0 


0.0 - 0.2 


10 - 100 


0.2 - 


0.5 


15 - 


60 




5kpc-II 


0.7 - 


1.0 


0.0 - 0.4 


< 50 


0.8 - 


1.5 


40 - 


70 


14060-1-2919 


5kpc 


0.5 - 


0.7 


0.4 


7 


0.4 - 


0.5 


60 - 100 




Ap A 


0.5 - 


2.0 


0.0 - 0.2 


< 100 


0.0 - 


0.9 


25 - 


55 




Ap B 


0.3 - 


0.7 


0.0 - 0.4 


7 - 30 


0.1 - 


0.4 


45 - 


90 




Ap C 


0.3 - 


0.7 


0.0 - 0.4 


4 - 7 


0.2 - 


0.4 


45 - 


82 


14348-1447 


5kpc-I 


0.3 - 


0.7 


0.0 - 0.6 


5 - 20 


0.4 - 


0.7 


50 - 


80 




5kpc-II 


0.3 - 


2.0 


0.0 - 0.4 


7 - 8 


0.4 - 


0.6 




40 












50 - 80 


0.0 - 


0.3 


50 - 


80 




Ap A 


0.3 - 


2.0 


0.0 - 0.4 


5 - 60 


0.0 - 


0.5 


45 - 100 




Ap B 


0.5 - 


1.0 


0.0 - 0.4 


< 30 


0.5 - 


1.2 


35 - 


60 




Ap C 


0.3 - 


2.0 


0.0 - 0.4 


< 100 


0.0 - 


0.9 


14 - 


80 




Ap D 


0.5 - 


0.7 


0.2 - 0.4 


4 - 50 


0.0 - 


0.6 


20 - 


80 


14394-1-5332 


5kpc-I 


0.5 - 


0.7 


0.2 


6 - 7 




0.7 




75 




5kpc-II 


2.0 




0.0 - 0.2 


6 - 100 


0.0 - 


0.8 


25 ■ 


-40 




A 


0.5 - 


2.0 


0.0 - 0.4 


4 - 100 


0.0 - 


0.5 


30 - 


75 




B 


0.7 - 


2.0 


0.0 - 0.2 


7 - 40 


0.3 - 


0.5 


30 - 


55 


15130-fl958 


5kpc 


0.3 - 


1.0 


0.0 - 0.4 


< 100 


0.0 - 


1.3 


20 - 


70 


15206+3342 


5kpc 


0.5 - 


0.7 


0.5 - 0.6 


1 - 3 


0.4 - 


0.5 


70 - 


80 






0.3 - 


2.0 


0.0 - 0.4 


10 




0.3 




91 



Table 9. Modelling results for the CS using Comb III (lYSP -I- VYSP). Col (1): object name. Col (2): aperture label as indicated in 
Figure 12] Col (3): range of ages for the lYSP component that adequately fit the data. Col (4): range of E{B — V) values for the lYSPs 
in Col (3). Col (5): range of ages of the VYSP that adequately fit the data. Col{6): range of E{B ~ V) values for the VYSPs in Col (5). 
Col (7): flux contributions of the VYSP components to the model, in the normalising bin. 
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Object 


Aperture 


Age of 


E(B • 


-V) 


Age of 


E(B 


-V) 


%VYSP 


Name 




lYSP 


of 


VYSP 




of 




IRAS 




(Gyr) 


lYSP 


(Myr) 


VYSP 




(1) 


(2) 


(3) 




(4) 


(5) 




(6) 


(7) 




5kpc 


0.5 - 


0.7 


0.4 


< 60 


0.8 - 


1.3 


16 - 45 




Ap A 


0.6 - 


0.7 


0.0 - 


0.2 


< 100 


0.0 - 


0.2 


< 5 




Ap B 


0.5 - 


0.7 


0.0 - 


0.2 


< 100 


0.0 - 


0.3 


< 20 




Ap C 


0.5 - 


0.7 


0.3 - 


0.5 


< 100 


0.0 - 


1.0 


5-15 




Ap D 


0.5 - 


0.9 


0.0 - 


1.0 


< 80 


0.9 - 


1.8 


30 - 63 




Ap E 


0.5 - 


0.9 


0.2 - 


0.6 


< 60 


0.9 - 


1.6 


22 - 47 




Ap F 


0.5 - 


0.7 


0.2 - 


0.6 


< 80 


0.0 - 


1.5 


5-36 




Ap G 


0.5 - 


0.7 


0.0 - 


0.3 


< 70 


0.0 - 


0.7 


5-20 




ApxOTAL 


0.5 - 


0.7 


0.2 - 


0.6 


< 100 


0.0 - 


1.5 


5-35 


15Ji27+234nCPA75") 


Ap A 


0.5 - 


0.7 


0.2 - 


0.5 


< 100 


0.0 - 


2.0 


< 20 




Ap B 


0.5 - 


0.7 


0.2 - 


0.6 


< 100 


0.0 - 


1.0 


5-25 




Ap C 


0.5 - 


0.7 


0.2 - 


0.6 


< 50 


0.0 - 


1.2 


5-30 




Ap D 


0.5 - 


0.7 


0.2 - 


0.6 


< 60 


0.0 - 


1.5 


5-25 




Ap E 


0.5 - 


0.7 


0.2 - 


0.5 


< 50 


0.0 - 


2.0 


< 30 




Ap F 


0.6 - 


0.9 


0.1 - 


0.4 


< 100 


0.0 - 


1.2 


< 25 




Ap G 


0.7- 


0.9 


0.0 - 


0.3 


< 100 


0.0 - 


0.8 


< 25 




ApTOTAL 


0.5 - 


0.7 


0.2 - 


0.5 


< 100 


0.0 - 


1.2 


5-25 


15Ji274-234nCPA75*'l 


Ap A 


0.5 - 


0.7 


0.0 - 


0.3 


< 60 


0.0 - 


0.6 


10 - 35 




Ap B 


0.5 - 


0.7 


0.0 - 


0.3 


< 100 


0.0 - 


1.5 


5-25 




Ap C 


0.5 - 


0.7 


0.0 - 


0.3 


< 100 


0.0 - 


1.5 


< 40 




Ap D 


0.5 - 


0.8 


0.0 - 


0.3 


< 60 


0.0 - 


1.5 


< 40 




Ap E 


0.6 - 


0.9 


0.0 - 


0.3 


< 60 


0.0 - 


1.0 


< 25 




Ap F 


0.6 - 


0.9 


0.0 - 


0.3 


< 100 


0.0 - 


2.0 


< 15 




ApxOTAL 


0.5 - 


0.8 


0.0 - 


0.3 


< 60 


0.0 - 


1.5 


5-30 


1 fi4744-'H4'Hn 




0.5 - 


0.7 


0.0 - 


0.4 


7-8 


0.3 - 


0.4 


65 - 80 




51{nr-TT 

tJlV^L/ J. J. 


0.5 - 


2.0 


0.0 - 


0.2 


6-8 


0.2 - 


0.4 


30 - 50 














50 - 90 


0.0 - 


0.2 


50 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


6-8 


0.3 - 


0.5 


40 - 75 














30 - 40 


0.3 - 


0.5 


40 - 75 


1 fi4874-'i447NF, 


5kpc 


2.0 




0.1 


< 80 


0.0 - 


0.8 


10 - 20 




Ap A 


2.0 




0.0 - 


0.2 


5-8 


0.2 - 


0.3 


45 - 55 














30 - 40 


0.0 - 


0.2 


50 - 60 


164874- 5447NW 


5kpc 


0.5 - 


0.7 


0.0 - 


0.4 


4-20 


0.1 - 


0.3 


55 - 85 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


4-60 


0.0 - 


0.4 


15 - 65 




Ap B 


0.7 - 


2.0 


0.0 - 


0.4 


4-8 


0.2 - 


0.3 


65 - 80 


1 7028+ 'i81 7NF, 




0.5 - 


1.0 


0.0 - 


0.4 


4-8 


0.3 - 


0.5 


50 - 65 


1 70284-581 7NW 


5kpc 


0.5 - 


1.0 


0.0 - 


0.2 


< go 


0.8 - 


1.5 


40 - 70 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


< 100 


0.0 - 


0.9 


30 - 85 


20414-1651 




0.5 - 


2.0 


0.0 - 


0.4 


4 - 100 


0.4 - 


0.8 


30 - 65 


21208-0519 


5kpc I 


0.5 - 


2.0 


0.0 - 


0.4 


5-50 


0.3 - 


0.5 


50 - 82 




5knr TT 


- 












- 






Ap A 


2.0 




0.0 


70 - 100 


0.3 - 


0.4 


40 - 50 




Ap B 


0.5 - 


2.0 


0.0 - 


0.4 


4-8 


0.0 - 


0.6 


35 - 70 














30 - 100 


0.0 - 


0.6 


35 - 70 


22491-1808 


5kpc 


0.3 - 


0.5 


0.0 - 


0.4 


4-8 


0.2 - 


0.4 


45 - 75 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


<100 


0.0 - 


0.8 


25 - 57 




Ap B 


0.3 - 


0.7 


0.0 - 


0.4 


4-7 


0.3 - 


0.4 


60 - 90 




Ap C 


0.3 




0.2 - 


0.4 


5 - 6 


0.2 - 


0.4 


44 - 60 


23233-f2817 


5kpc 




















Ap A 




















Ap B 


















23234-1-0946 


5kpc I 


0.5 - 


2.0 


0.0 - 


0.4 


6 - 30 


0.4 - 


0.6 


60 - 82 




5kpc II 


1.0 - 


2.0 


0.0 - 


0.2 


< 100 


0.5 - 


1.5 


17 - 37 




Ap A 


0.7 




0.0 


< 20 


0.5 - 


0.8 


15 - 25 


23327-1-2913 


5kpc I 



















5kpc II 
Ap A 
Ap B 



Table 9. Continued 



The properties of the stellar populations in ULIRGs I: sample, data and spectral synthesis modelling. 
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Object 


Aperture 


Age of 


E{B ■ 


- V) 


Age of 


E(B -V) of 


%VYSP 


Name 




lYSP 


of 


VYSP 




of 






IRAS 




(Gyr) 


lYSP 


(Myr) 


VYSP 






(1) 


(2) 


(3) 




(4) 


(5) 




(6) 




(7) 


08572+3915 


5kpc 


0.5 - 


2.0 


0.0 - 


0.4 


6 - 70 


0.1 - 


0.6 


30 


- 90 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


< 100 


0.0 - 


1.0 


10 


- 50 




Ap B 


0.3 - 


2.0 


0.0 - 


0.2 


< 100 


0.0 - 


0.5 


10 


- 60 




Ap C 


0.3 - 


2.0 


0.0 - 


0.4 


< 100 


0.0 - 


0.7 


25 - 


100 


10190+1322 


5kpc I 


0.7 - 


2.0 


0.0 - 


0.4 


< 60 


0.7 - 


1.1 


30 


- 65 




5kpc II 


0.7 




0.2 - 


0.4 


50 - 100 


0.8 - 


1.1 


30 


- 55 




Ap A 


0.5 - 


0.7 


0.0 - 


0.2 


< 100 


0.0 - 


2.0 


< 25 




Ap B 


0.5 - 


1.0 


0.0 - 


0.4 


8 - 30 


0.8 - 


0.9 


52 


- 66 


10494+4424 


5kpc 


0.5 




0.0 - 


0.2 


5 - 10 


1.1 - 


1.3 


30 


- 45 














40 - 60 


0.8 - 


0.9 


50 


- 70 




Ap A 


- 








- 




- 








Ap B 


0.3 - 


0.7 


0.0 - 


0.4 


6 - 60 


0.0 - 


0.6 


45 


- 75 


13305-1739 


5kpc 


0.5 - 


2.0 


0.0 - 


0.4 


4 - 7 


0.4 - 


0.7 


50 


- 80 


14252-1550 


5kpc 


0.5 - 


2.0 


0.0 - 


0.4 


< 90 


0.5 - 


1.3 


35 


- 80 




Ap A 


0.5 - 


1.0 


0.0 - 


0.4 


< 100 


0.0 - 


2.0 


5 


- 80 


16156+0146 


5kpc I 


0.7 - 


2.0 


0.0 - 


0.4 


9 - 10 




0.5 


60 


- 80 




5kpc II 


0.3 - 


2.0 


0.0 - 


0.4 


< 100 


0.0 - 


0.6 


20 


- 70 




Ap A 


0.5 - 


2.0 


0.0 - 


0.4 


< 30 


0.0 - 


0.7 


15 


- 81 


17044+6720 


5kpc 


0.3 - 


1.0 


0.0 - 


0.4 


7 - 40 


0.4 - 


0.6 


40 


- 87 




Ap A 


0.3 - 


2.0 


0.0 - 


0.4 


< 100 


0.0 - 


0.8 


30 


- 85 




Ap B 


0.3 - 


2.0 


0.0 - 


0.4 


4 - 100 


0.1 - 


0.4 


40 


- 85 


17179+5444 


5kpc 


0.7 - 


2.0 


0.0 - 


0.4 


4 - 100 


0.3 - 


1.1 


25 


- 60 


23060+0505 


5kpc 






















Ap A 


1.0 - 


2.0 


0.0 


6 - 7 




0.7 


50 


- 60 














40 - 60 


0.4 - 


0.5 


50 


- 60 




Ap B 


0.7 - 


2.0 


0.0 - 


0.4 


10 - 50 


0.0 - 


0.5 


30 


- 75 


23389+0303 


5kpc I 


0.7 - 


1.0 


0.0 - 


0.4 


6 - 8 


0-6 0.8 


50 


- 72 




5kpc II 


0.7 - 


2.0 


0.0 - 


0.2 


< 100 


0.2 - 


0.9 


25 


- 60 




Ap A 


0.5 




0.0 


10 - 20 




0.7 




- 40 




Ap B 


0.7 - 


2.0 


0.0 - 


0.4 


7 




0.8 




75 




Ap C 


0.7 - 


1.0 


0.0 - 


0.4 


5 - 60 


0.2 - 


0.8 


35 


- 60 




Ap D 


0.7 - 


2.0 


0.0 - 


0.2 


4 - 6 


0.6 - 


0.7 


40 


- 50 














40 - 100 


0.2 - 


0.4 


45 


- 70 



Table 10. Same as Tabled but for the 10 additional ULIRGs included in the ES. 
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APPENDIX 
NOTES ON INDIVIDUAL SOURCES 

IRAS 00091-0738: optical and near-IR images reveal 
the presence of several bright knots sho wing the loca t ion o f 
young ( < 7 Myr) stellar populations (jSurace et all |2000| ) 
in this double nucleus system (NS = 2.2 kpc). The double 
nucleus structure is not resolved in our spectra and a single 
5 kpc aperture was used during the analysis. Our modelling 
results reveal the presence of a significant contribution from 
VYSPs in the nuclear region, which is consistent with the re- 
sults obtained for the extended apertures. IRAS 00091-0738 
is classified as an Hll-galaxy at optical wavelengths. There- 
fore, it is possible to use the nebular emission lines in order 
to provide a further constraint on the ages of the VYSP 
component. The VYSP age obt ained using the equiva- 
lent width and the iLe ithe rer et al. (1999) models (EW(Hq) 
= 28 A, tvYSP ~ 7 Myr) is consistent with the lower end 
of the range of ages found from the modelling. 

IRAS 00188-0856: this compact, single nucleus sys- 
tem is classi fied as a LINER at op tical wavelengths. On the 
other hand, ilmanishi et al.1 (|2007h and IValdes et all (120051 ) 
reported evidence of the presence of an AGN in the nucleus 
of the galaxy based on their near- and mid-IR studies of 
the object. However, the AGN does not dominate the IR 
emissi on, which is likely to be domin ated by starburst ac- 
tivity (|Lutz et al.l . ll999i : IValdes et al.L.2005. ). The modelling 
results suggest the presence of a dominant lYSP with a rel- 
atively small percentage contribution of any VYSP. Note 
that, due to the small contribution of the VYSP to the op- 
tical light from the galaxy, both the age and reddening of 
that component are relatively unconstrained. 

IRAS 01004-2237: the optical spectrum of this single 
nucleus system shows a mixture of HIT and Seyfert 2 features 
(|Veilleux et al.. 1999al). with clear evide nce of a buried AGN 
at IR wavelegths ( Imanishi et al.l . |2007| ). The modelling re- 
sults clearly suggest the presence of a dominant VYSP with 
relatively low reddening. The VYSP ages obtained from the 
modelling (tvYSP = 9-20 Myr) differs than those obtained 
from the Ha equivalent width (EW(Ha) = 58 A, tvYSP '~ 
6.5 Myr ). However, given the uncertainties attached to the 
synthetic templates for very young ages (< 10 Myr), we con- 
sider that ~ 3 Myr is not a major disagreement between the 
two results. 

IRAS 08572-1-3915: this is double nucleus system 
consisting of two spirals galaxie s in the process of merg- 
ing, rotating in different planes (|Arribas et al.l . |2000| ). The 
nuclear, 5 kpc aperture was centred in the northern nu- 
cleus (the secondary nucleus was not covered by the slit), 
t hough t to be responsible for the far-IR emission from the 
source (ICondon et al.l.[l99ll ; ISoifer et all l2000l : lEvans et all 
200^;'NaKar et al.V2003^. The /fS'r-WFPC2 image taken by 
Sur acc ct al. (1998) reveals the presence of several young 
{^knota = 5-10 Myr), massive knots located in both the 
tails and the nuclear regions. The optical spectrum of IRAS 
08572-1-3915 is a composite of HIT and LINER-like fea- 
tures. Although no evidence for the presence of an AGN 
is found in the optica l and near- IR (jVeilleux et al.. 1999al : 
IVeilleux et al.. 1999bl ). IRAS 08572-F3 915 shows signs of 
AGN activity at mid-IR wavelengths (jlmanishi fc Dudlevl 
l2002l : Ilmanishi etall . l2006l : lArmus et al.l . l200'iD . Moderately 
reddened VYSPs are found in the norhern nucleus, which 



contribute significantly to the optical light. The reddenings 
and ages of the stellar components in the extended apertures 
are consistent with those of the nuclear aperture. 

IRAS 10190-1-1322: The north-east galaxy of this 
double nuleu s system is b ri ghter at near-IR and radio 
wavelengths (iMurphv et al.1 . 1200 ll : iDasvra et ail . l2006bl : 
iNagar et al. 1. 120031 ) ■ while the south-wes tern nucleus appears 
the brightest in the optical. Recently, iRupke et al.l (|2002l ) 
classified the north-east nucleus as a LINER and the south- 
western source as an HII galaxy in their high-resolution spec- 
troscopic study of ULIRGs. No evidence for AGN activity 
is found at near- and mid-IR wavelengths (Veilleux et al., 
1999b, Imanishi et al., 2007), which is consistent with the 
idea of starburst activity powering the IR luminosity of the 
source (Farrah et al., 2003). The results obtained from the 
modelling suggest that the VYSPs are concentrated in the 
two nuclei of the galaxy, and in the region in between, coin- 
ciding with large concentrations of gas and dust. In this case, 
only the extracted spectrum from aperture 5kpc-I shows 
clear Hll-like region features. The VYSP age derived us- 
ing the Ha equivalent width (EW(Ha) — 47 A, twsp ~ 
6.5 Myr) for that region of the galaxy is consistent with the 
lower end of the range of ages found for the VYSPs from 
modelling the optical continuum. 

IRAS 10494-1-4424: the modelling resuhs obtained 
using near-IR and optical photometric points clearly sug- 
gest a s tarburst as an ioniz ing source for this single nucleus 
system (iFarrah et al.l. |2003|). which is classified as a LINER 
in the optical ( Kim &: Sandersl . ll998l ). Modelling the optical 
spectra, we find that the VYSPs located in the nuclear and 
the northern region of the galaxy make a significant contri- 
bution to the optical light. However, adequate fits for Ap 
A (sampling a region towards the south of the galaxy) are 
only obtained using Comb I, with the optical light domi- 
nated by a 400 Myr lYSP with little or no reddening, plus a 
small contribution from a 12.5 Gyr OSP. In terms of redden- 
ing, our modelling results clearly provide evidence for higher 
reddenings towards the nuclear regions of the galaxy. 

IRAS 12072-0444: this ULIRG appears as a single nu- 
cleus system in the optical and near-IR images of iKim et al.l 
(2002). However a second nucle us becomes apparent in the 
i/ST-NICMOS H-band images (|Veilluex et al.l . |2006| ). This 
second nucleus was first detected in the WFPC2 images of 
ISurace et al. I 11998), bu t was cla ssified as a super star clus- 
ter. More recently, Dasvra et al.l ( j2006a ) confirmed spectro- 
scopically the presence of a second genuine nucleus, with 
a NS = 2.8 kpc. The double nucleus structure is unre- 
solved in our spectroscopic data. The detection of strong 
[Si IV] emission is evidence for the presence of a buried 
AGN in this source , class i fied as a Sy2 ga l axy i n the op- 
tical (IVeilleux et al ], 1 19951 : iKim fc Sand"er3. 1 19981) . Moder- 
ately reddened VYSPs which contribute significantly to the 
optical light from the galaxy are located in the nuclear re- 
gion. Comparing the nuclear with the extended apertures, 
the modelling results suggest the presence of younger, more 
highly reddened VYSPs in the nuclear regions of the galaxy 
than in the extended apertures. 

IRAS 12112-1-0305: this double nucleus system is 
classifies as a LINER in the optical, shows no evidence for 
powerful AGN activity (jCenzel et al.|[l998l : Veilleux et al., 
1999b), and is classified as a starburst at all other wave- 
lengths studied (|Lutz et al.l . Il999l : iFranceschini et al] . l2003l : 
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iRisaliti et al] , l2006l : lArmus et all. l2007h . The modelling re- 
sults show the presence of dominant, relatively high red- 
dened VYSPs located in the north-eastern nucleus (5 kpc- 
I), which is the brightest at optical wavelengths. For this 
object, the ages and percentage contributions to the flux for 
both nuclei and the region in between are consistent. On the 
other hand, our results suggest increasing concentrations of 
obscuring dust towards the northern source. 

IRAS 12540+5708: also known as Mrk231, IRAS 
12540+5708 is the most luminous ULIRG in our sample 
(log(LiR) = 12.54). This object is classified as a Syl galaxy 
in the optical, and considered an AGN -powered ULIRG 
based on mid- and near -IR spectroscopy (|Lutz et al.l . 1 19991 : 
iRigopoulou et al.l . ll999l ). However, the fraction of the bolo- 
metric luminosity produced by the AGN remai ns controver- 
sial. D espite its warm mid- to far-IR colours, iDavies et al.l 
l|2004D found that star formation within the inner 300pc 
is responsible for 25% - 40% of the entire luminosity of 
the galaxy, and that starburst activity seems to power 
most of the far-IR luminosity l|Downes fc Solomonl . 1 19981 : 
iFarrah et all |2003| ). On the other hand, ages from a few 
Myr in the inner parts of the galaxy, up to ages of or- 
der of Gyr in the extended regions, have been found for 
the stellar populations in the galaxv ([Surace et all . 1 19981 : 
ICanalizo fc Stocktonl l2000al : iDavies et al.1 . l20oi ) 

Due to ghost images affecting many of the apertures 
in the red arm, it was only possible to model the extracted 
spectra from the blue arm, with the exception of Ap C. In 
addition, the powerful AGN emission prevents any attempt 
at modelling the stellar populations for those apertures sam- 
pling the nuclear regions (5 kpc and Ap D). Therefore, only 
the extracted spectra from apertures A, B, C and E were 
used for the analysis of the stellar populations presented 
here. Note t hat Ap C s amples the bright condensation (the 
"horseshoe" . iHamilton fc Keel, 1987) ^ 3 kpc to the south 
of the galaxy , as se en in the HST-WFPC2 images taken by 
ISurace et al.l (|l998l '). Figure [7] shows the positi on of the slit 
and th e extent of Ap C, superimposed on the lSurace et al.l 
l|l99 ^ _ffS'r-WFPC2 F439W image. In order to minimize 
the potential AGN contamination (including the seeing disk 
of the bright AGN and any dust-scattered AGN light), the 
spatial axes of the blue and red 2-D frame were reversed. 
The reversed and the original frames were registered and 
then subtracted. A spatial cut of the subtracted 2-D frame 
is shown in Figure [8l as well as the spatial cut of the origi- 
nal frame and the extent of Ap C. It is clear from the figure 
that any AGN contamination has been removed, and there- 
fore, we were able to model the full spectral range for this 
aperture. 

The modelling results for Ap C reveal the presence of a 
VYSP with a low reddening {E{B-V) = 0.0 - 0.1) and ages 
in the range 50 < ivYSP < 80 Myr, that accounts for almost 
the entire optical light within the aperture (Pvysp = 70 - 
96%). It is worth noting that Comb I, Comb II and Comb 
III provide consistent results for this aperture. In addition, 
the extracted spectrum shows the clear presence of He I ab- 
sorption features at various wavelengths, which are indicated 
in Figure |9l In the models, the He I lines are strongest for 
single stellar populations (SSP) with ages in the range 20 
- 50 Myr, when the stellar populat ion is dominated by B 
stars l|Gonzalez Delgado et al.l . l2005l ). and are not observed 
in stellar populations older than 100 Myr (the lifetime of 
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Figure 7. Mrk 231: the slit position and t he location of Ap C 
are shown in the figure superimposed on the ISurace et al.l (|l998l ') 
F439W WFPC2 HST image. The major ticks are arcseconds. 



IRAS 12540+5708 (Mrk231) 




Figure 8. Mrk 231: spatial cuts (4400 - 4600 A) extracted from 
the subtracted and original 2-D frames showing the region includ- 
ing the bright condensation to the South of the galaxy. It is clear 
from the figure that the AGN contamination has been almost en- 
tirely removed after applying the technique described in the text. 
The location of Ap C is also shown in the figure. 



B stars). Therefore, the presence of such features helps to 
provide a further constraints on the age of the VYSP com- 
ponent within the aperture. The detection of the He I ab- 
sorption lines is consistent with the lower end of the range 
of ages of the VYSP component deduced from modelling the 
continuum. 

Previous spectroscopic studies already revealed the 
presence of YSP s in th e host galaxy of M rk 23 1 
(Hamilto n fc Keell, Il987l: ICanalizo fc Stocktonl l2000bl ). 
[Canalizo fc Stocktonl ( 2000b ) carried out a similar study to 
that presented here, modelling the optical spectra extracted 
from a series of apertures sampling various regions of the 
galaxy. In order to fit the extracted spectra, they used a 
combination of a 10 Gyr OSP plus unreddened YSPs with 
varying ages. This is similar to the Comb I used here, al- 
though reddening effects are only fully considered in this 
study. Their "regions" a, b, c and k are similar to our Ap C, 
Ap B, Ap A and Ap E, respectively, shown in Figure[51 They 



found YSPs with ages of 5 Myr, 140l^^ Myr and 42 Myr 
for region c, b and a, respectively, which is consistent with 
the modelling results found here using Comb I. On the other 
hand, they found no evidence for significant YSPs in region 
k, similar to the Ap E used here. However, we find impor- 
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Figure 9. Extracted spectrum from Ap C for Mrk 231. The figure 
shows the detection of the He I absorption lines at 3820 A, 4026 A, 
4388 A and 4471 A. These hues have a maximum strength for 
ages in the range of 20 — 50 Myr, and are not produced by stellar 
populations older than 100 Myr. 



tant YSPs (YSP% = 50 - 70%) within that aperture, with 
ages and reddening in the range of tvsp = 60 - 100 Myr and 
E{B -V) = 0.2 - 0.5. This difference is likely to be due to 
the fact that the rest wavelength range of the spectra used 
in their study is 3800 - 7000 A. Therefore, they do not sam- 
ple the near-UV, which is a crucial region of the spectrum 
for constraining the properties of the YSPs. In addition, it is 
notabl e that only unredden e d synth etic templates are used 
in the lCanalizo fc StocktonI ()2000bl ) study, whilst reddened 
YSPs are required to fit the data in the case of Ap E. 

IRAS 13305-1739: the optical spectrum of this very 
compact ULIRG is that of a Sy2 galaxy. In addition, 
IVeilleux et al.. 1999bl found broad Paa and [Si IV] emission 
at near-IR wavelengths - clear evidence of the presence of a 
quasar/AGN - in the nucleus of the galaxy. A single nuclear 
aperture was used for the analysis of the stellar populations 
in this object. We find a moderately reddened, dominant 
VYSP component with a remarkably narrow range of ages 
(4 < tvYSP < 7 Myr). 

IRAS 13428+5608: also known as Mrk 273, this 
galaxy shows an impressive tidal tail extended over 30 kpc 
towards the south of the gal axy. Seen as a single - nucleus ob- 
ject in UV/optical images (jSurace et al.l , I2OO0I : iKim et al] . 
I2OO2I'). a double nucleus structure emer ges at near-IR 



, Armus et al. ^ 1990l: iMaiewski et al.l . Il99 Knapen et al, 
1997; Scoville et al.l . bOOd) and rad io (lUlvestadl 1 1984 
"Condon^eTaLl . Il99ll : ICole et all . Il999l l wavelengths, with a 
nuclear separation of ~ 700 pc. The double nucleus struc- 
ture is unresolved in our spectroscopic data. Classified as 
a Sy2 at optical wavelengths, IRAS 13428-f5608 is a com- 
posite object power ed by both starburst and AGN activity 
j Armus et al.l ■ l2007^ . However, the dominant power source 
remains controversial. 

In order to perform a detailed study of the stellar pop- 
ulations detected at optical wavelengths, a total of 7 aper- 



tures were extracted, sampling the nuclear region and the 
large tidal tail. In the case of the nuclear 5 kpc aperture, we 
find an "old" lYSPs (Iiysp = 0.7 - 2.0 Gyr) with low red- 
dening which dominates the optical emission, with a smaller 
contribution from a VYSP component with higher redden- 
ing. Previous optical spectroscopy already revealed the pres- 
ence of a mix of stellar populations in t he nuclear region of 
Mrk273 (jGonz alcz Dclg ado et ahl . l200lh . In a broad-brush 
se nse, the results presented here are consistent with those 
of iGonzalez Delgado et al] (I2OOII ). The extended apertures 
sample the prodigious tidal tail to the south of the galaxy 
and the north of the nuclear region. Overall, the results 
show a remarkable uniformity, with a dominant "old" lYSP 
present at all locations in the galaxy, out to a radius of 30 
kpc. It is also notable that a VYSP component is required in 
order to fit the data for all the apertures, even for those sam- 
pling the extended tail. No evidence suggesting the presence 
of important age or reddening variations across the body of 
the galaxy is found in the case of Mrk 273. 

IRAS 13451+1232: also known as PKS1345+12, this 
is the only ULIRG in our sample also classified as a radio 
galaxy. The activity is concentrated in the western nucleus, 
which is the only nucleus covered by the two slit positions 
used for this o bject. This nucl eus is classified as a Sy2 in the 
optical ( V eilleux et al.. 199931), and a compact source at ra- 
dio wavelengths ( Evans et al.l , 1 19991 ) . Significant young stel- 
lar populations have already been dete cted in PKS1345+12 
at optical and near-I R wavelengths jSurace et al.l . Il998l : 
iTadhunter et al.l |2005| ) , which may have been formed as a 
consequence of the merger. 

A detailed analysis of the stellar populations within 
the galaxy, using the s o-called Comb I, was a l ready pre- 
sented for this object in [Rodriguez Zaurm et al] (|2007l ). We 
emphasize that prior to mod elling the data for both the 
iRodrfguez Zaurm et al.l (|2007l ) study and this work, we sub- 
tracted the best model for the nebular continuum emission 
from the detailed emission line modelling work of Holt et al. 
(2003), in the case of the 5 kpc aperture (i.e. the "NUC" 
aperture in Rodriguez Zaurfn et al., 2007). Such emission 
is also likely to be important in the cases of Ap A and Ap 
B in PKS 1345+12. For these two apertures, both the re- 
sults from the modelling using the nebular corrected and 
uncorrected spectra are presented in Table [5] 

In this study, we use an updated version of the code 
that allows a more detailed analysis of the different, impor- 
tant regions of the optical spectra and, therefore , it is possi- 
ble to refine the results presented in [Rodriguez Zaurin et al.l 
(|2007l ). In addition, the extracted spectra are also modelled 
using Comb II and HI. In general terms, the results found 
here are consistents with those of [Rodriguez Zaurm et al.l 
iiool), suggesting that a VYSP may not be as important 
in the case of PKS 1345+12 as for other ULIRGs. In ad- 
dition, we find that adequate fits are also obtained for all 
the apertures with an "old" lYSPs (tivsp = 1.0 - 2.0 Gyr) 
that dominate the flux at all locations sampled by the aper- 
tures. Due to the small percentage contribution of the VYSP 
component, the age and reddening of this component are rel- 
atively unconstrained, which explains the wide range of ages 
and reddening presented in Table [5] Finally, note that this 
is the most massive galaxy in our sample, in terms of its 
stellar population. 
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IRAS 13539+2920: two 5kpc apertures, sampling the 
two nuclei of this double nucleus system, were extracted for 
this object. For the western nucleus (5kpc-II), which is the 
brightest at optical wavelengths, we find highly reddened 
VYSPs that contribute significantly or dominate the optical 
emission. For the eastern nucleus (5kpc-I), we find ages con- 
sistent with those of the western source, but with lower red- 
dening values. On the other hand, IRAS 13539-1-2920 is spec- 
troscopically classifed as an Hll-galaxy in the optical. How- 
ever, modelling the emission lines of the extracted spectra 
for aperture 5kpc-II, sampling the nucleus to the west of the 
galaxy, we find [OIII]5007/H/3 and [NII]6583/Ha emission- 
line ratios consistent with a combination of an Hll-galaxy 
and a LINER. Moreover, the ID-spectrum of that nucleus 
shows clear signs of line splitting, indicating the presence 
of ionized gas outflows. On the other hand, no evidence of 
Hll-like region features is found in the case of the eastern 
source, sampled by aperture 5kpc-I (it has a LINER spec- 
trum). Therefore, no attempt to estimate the age of the 
VYSP using the Ha equivalent width was made for this 
galaxy. 

IRAS 14060+2919: we find a moderately reddened, 7 
Myr VYSP component that dominates the optical emission 
in the nuclear region of this single nucleus system. In general 
terms, the results for the extended apertures are consistent 
with those in the nuclear region, except for perhaps a slightly 
smaller VYSP contribution in Ap A (although Ap A is not 
so well constrained because it is noisier). In addition, IRAS 
14060+2919 is classifed as an HII galaxy at optical wave- 
lengths. The Ha equivalent widths measured for Ap B, Ap 
C and 5kpc ( 77 A < EW{Ha) < 177 A) reveal the presence 
of VYSPs younger than ~ 6 Myr, which is consistent with 
the continuum modelling results. 

IRAS 14252-1550: classified as a LINER in the op- 
tical, no signatures of buried quasar activity are found in 
the 5-35 nm mid-IR s pectrum of this double nucleus system 
l|lmanishi et al.l . |2007| ) . The modelling results are consistent 
with a wide range of VYSPs ages, reddenings and percent- 
age contributions. No clear evidence suggesting the presence 
of important age or reddening variations is found. 

IRAS 14348-1447: This is a double nucleus system, 
in which the ULIRG activity is likely associated with the 
southern nucleus, which dominates the near-IR luminos- 
ity. Several bright knots are seen in t he in the U'-, B- 
and I-band images l|Surace et al urace fc Sander j . 

l200d ). some o f the knots with e stimated ages younger 
than 5-7 Myr (|S urace et al.l . l200Cll V The o ptical spectrum 
of IRAS 14348-1447 is th at of a LINER l|Veilleux et al.l . 
1 19951 : iKim fc Sanderg, Il998l ) , and is classified as as starburst- 
powered ULIRG at n ear-I R, mid-IR , and X-ra y wavel e ngths 
llVeilleux et al.l.[l9 97: Gcnz el et ai].il9 98: Lut z et al.l.ll999l : 
Franceschini et al.1 . 120031 : iFarrah et al.1 . 120031 : IRisaliti et all . 



20061 ). We find significant, or even dominant VYSPs, with 



a wide range of reddening values, located in the nuclei of 
the system. Note that in the case of the north-eastern nu- 
cleus (5kpc-II) the model solutions for the VYSPs divide 
into two groups: a younger group of 7 - 8 Myr, and a sec- 
ond, older, group with ages of 50 - 80 Myr. In this case, 
we cannot distinguish between the two groups using the de- 
tailed fits. However, modelling the emission lines of the ex- 
tracted spectra for aperture 5kpc-II we find [OIII]5007/H/3 
and [NII]6583/Ha emission-line ratios consistent with an 



Hll-galaxy. The VYSP age obtained using the Ha equiva- 
lent width (EW(Ha) = 51 A, tvYSP ~ 6.5 Myr) is consistent 
with the first group of adequate fits, including the younger 
YSP component. Finally, no evidence is found for significant 
variations in the YSP age and reddening across the body of 
the galaxy. 

IRAS 14394-1447: this IRAS source is a muhiple (> 
2 nucleus) system, with the two main components separated 
by 54 kpc - the widest nuclear separation in the ES. The 
eastern nucleus is classified as a Sy2 and identified with 
the source responsible for the far-IR luminosity (Veilleux et 
al., 1999a, Rupke et al., 2005). Dominant, highly reddened 
VYSPs are found for that nucleus (5kpc-I). On the other 
hand, the modelling results are consistent with a wide ranges 
of age and reddening for the VYSP component in the case of 
the 5 kpc aperture extracted for the western source (5kpc- 
II). Moreover, adequate fits are only obtained for that aper- 
ture using a 2 Gyr lYSP which, in fact, dominates the fiux 
contribution in the normalising bin. These results clearly 
demonstrate large differences between the stellar popula- 
tions in the two nuclei, perhaps indicating the presence of a 
more powerful starburst activity in the eastern component. 
In addition, two further apertures were extracted to the east 
and west of the western source. We find no clear evidence for 
the presence of large age and/or reddening variations across 
the body of that source. 

IRAS 15130-1958: the modelling results are consis- 
tent with VYSPs with a wide range of ages, reddenings and 
percentage contributions for the single 5kpc aperture ex- 
tracted for this single nucleus system, classified as a Sy2 
in the optical. In addition, Wolf-Rayet (WR) features are 
detected for this galaxy. The detection of such features in- 
dicates ongoing star formation activity. Figure [TOl shows the 
wavelength range 4450 - 4750 A for the spectrum of the 
galaxy and a mid-type WN5-6 star spectrum in the Milky 
Way (Paul Crowther, private communication). The aim of 
this figure is to compare the morphology of the two spec- 
tra in the relevant wavelength range. The figure shows the 
detection of the NIII4634 and the HeII4686 features in our 
spectrum, consistent with the presence of WR stars. The de- 
tection of such features indicates ongoing star formation ac- 
tivity in the case of IRAS 15130-1958. Among all the galax- 
ies in our sample, this is the only case with a clear detec- 
tion of WR features. Assuming that the VYSP component 
is associated with the WR features, the detection of these 
features provides a further constraint on the VYSP ages, 
i.e. the VYSP ages must be < 10 Myr. This is consistent 
with the lower end of the range of VYSP ages obtained from 
modelling the continuum. 

IRAS 15206-3342: bright knots are observed in the 
innermost region of this single nucleus system, and also 
spread over even larger scales, spanning a wide range of ages 
(from a few Myr up to 2 Gyr, Suracc et al., 1998). This ob- 
ject is spectroscopically classified as an Hll-like galaxy in 
the optical (Kim et al., 1998; Veilleux et al., 1999a). More- 
over, no evidence for the presence of a buried AGN /quasar 
is found at I R and radio wavelengths (I Veilleux et al.l . Tl997l : 
iNagar et all [2003: Imanishi_ct_al., 2003) ■ The large nuclear 
Ha equivalent width and luminosity are a lso consistent with 
a sta rburst-dominated ionizing source (|Arribas fc Coliiial . 
l2002h . We find adequate fits for two groups of combinations 
of stellar populations. The VYSPs in the first group are 
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Figure 10. Comparison between the spectrum of IRAS 15130- 
1958 (top) and a mid-type WN5-6 star spectrum in the Milky 
Way (bottom). The figure shows the detection of the NIII4634 
and the HeII4686 features, revealing the presence of WR stars in 
the galaxy. 



younger than 3 Myr, while such component has a age of 
10 Myr in the case of the second group. However, the large 
equivalent width of the Ha emission line (EW(Hq) = 315A), 
suggest VYSP ages of tvYSP ^ 5 Myr. This result favours 
the first group of adequate fits, including the younger VYSP 
component, whic h is also consi s tent w ith the upper limit of 
20 Myr found by iFarrah et al] l|2005l ) on the basis of their 
UV/optical spectroscopic study of the stellar populations in 
the galaxy. 

IRAS 15327+2340: also known as Arp 220, this dou- 
ble nucleus system is the closest ULIRG in our sample 
(z = 0.018). The double nucleus structure is unresoved in 
our spectroscopic data. A biography and a detailed anal- 
ysis of modelling resul t s are presented for this object in 
iRodn'guez Zaun'n et al.1 (|2008'). For completeness, the mod- 
elling results of Comb I, Comb II and Comb III are shown 
in Table [5l Table [7] and Table |9l respectively. The extrac- 
tion apertures and the extracted spectra are also shown in 
Figures [5] and [5] 

IRAS 15462-0450: due to the powerful AGN emis- 
sion, it was not possible to model the stellar populations 
in this compact ULIRG, which is classified as a Syl galaxy 
at optical wavelengths. The extraction aperture and the ex- 
tracted ID-spectrum are shown in Figures [2] and |6] 

IRAS 16156+0146: three apertures were extracted 
for this double nucleus system, sampling the two nuclei and 
the region in between. The optical spectrum is that of a Sy2 
galaxy, with the AGN activity associated with the north- 
western source. Moderately reddened, dominant VYSPs of 
age 9-10 Myr are found for that source (5 kpc-I). On the 
other hand, in the case of the south-eastern nucleus (5 kpc- 
II), the modelling results are consistent with wide ranges 
of age, reddening and percentage contribution for both the 
VYSPs and the lYSPs. It is notable that the minimum per- 
centage contribution of the VYSPs found for the aperture 
5kpc-I (60%), is significantly higher than in the case of Ap 
A and 5kpc-II (20% and 15%, respectively). These results 
suggest higher concentrations of VYSPs towards the north- 



west of the system, although we cannot rule out the idea of 
significant VYSPs at all locations. 

IRAS 16474+3430: three apertures were used for the 
analysis of the stellar populations in this double nucleus sys- 
tem, centered on the nuclei and the region in between. The 
modelling results found for the southern nucleus (5kpc-I), 
the brightest at optical wavelengths, reveal the presence of 
moderately reddened VYSPs (tvYSP = 7-8 Myr) that dom- 
inate the flux contribution at optical wavelengths. For the 
other two apertures, the VYSPs are divided into two groups: 
a younger group with ages in the range 6-8 Myr, and a 
second, older, group with ages in the range 50 - 90 Myr. 
However, IRAS 16474+3430 is classified as an Hll-galaxy, 
exhibiting Hll-like region features at all locations covered 
by the slit. The values of the Ha equivalent widths ( 101 A, 
49 A and 69 A for 5kpc-I, 5kpc-II and Ap A, respectively) 
are only attained for stellar populations of ~ 6.5 Myr. This 
is consistent with modelling results found for the southern 
nucleus (5kpc-I), and the younger group of VYSPs in the 
cases of 5kpc-II and Ap A. 

IRAS 16 487+5447: this galaxy is classified as a 
LINER galaxy l|Kim fc Sanders! . [l998, ) in the optical, and 
shows no evidence for a buried AGN at mid-IR wavelengths 
(|lmanishi et al. 



200'iD. consistent with a starburst as the 



power source ( Lutz et al.l . Il999l ). Two slit positions were 



used during the obervations of this double nucleus system. 
For the eastern nucleus, the modelling results are consis- 
tent with the combination of a dominant "old" lYSP plus a 
VYSP with a wide range of ages and reddenings. It is no- 
table that adequate fits are also obtained for this nucleus 
using Comb I, with a dominant 12.5 Gyr OSP. On the other 
hand, for the western source, we find VYSPs with relatively 
low reddening, which represent a large fraction of the op- 
tical emission from the galaxy. In addition, modelling the 
emission lines of the 5 kpc aperture for this source we find 
[OIII]5007/H/3 and [NII]6583/Ha emission-line ratios con- 
sistent with an Hll-galaxy. The VYSP age obtained using 
the Hq equivalent width (EW(Ha) = 28 A, Ivysp ~ 7 Myr) 
is consistent with the lower end of the range of ages found 
from the modelling. These results clearly demonstrate large 
differences between the stellar populations in the two nuclei. 
Finally, no clear evidence is found for the presence of age or 
reddening trends in the case of the western source. 

IRAS 17028+5817: again, two slit PAs were used 
during the observations of this system, sampling the two, 
widely separated (25 kpc) galaxies. We find that VYSPs 
make a large contribution, or even dominate, the flux at op- 
tical wavelengths in the two nuclei. However, in the case of 
the eastern galaxy, classified as an Hll-galaxy, we find mod- 
erately reddened VYSPs that span a narrow range of ages. 
On the other hand, the modelling results are consistent with 
a wider range of ages and significantly higher reddening val- 
ues for the VYSPs located in the western galaxy, which is 
classifed as a LINER in the optical. The modelling results 
obtained for the eastern source are consistent with the age of 
~6 Myr estimated using the Ha equivalent width (EW(Ha) 
= 102 A). In addition, we extracted a third, extended aper- 
ture (labelled as Ap A in Figure [2}, sampling the region to- 
wards the south of that source. We do not find clear evidence 
suggesting the presence of large age or reddening variations. 

IRAS 17044+6720: classified as a LINER in the opti- 
cal, this single nucleus system shows strong evidence for the 
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J iresence of obs c ured AGN activity at mid-IR wavelengtlis 
Imanisfii et alJ . |2007| ) . Tlie modelling results found in the 
nuclear region of the galaxy are consistent with VYSPs of 
relatively low reddening, that make a large contribution to 
the fliix at optical wavelengths. Comparing these results 
with those found for the extended apertures, we do not find 
clear evidence for the presence of significant age and/or red- 
dening variations across the body of the galaxy. Note that 
Ap B (see Figure O has line ratios typical of an HII re- 
gion ([NII]6583/Hq = 0.24 and [OIII]5007/H/3 = 3.0). The 
estimated Ha equivalent width is EW(Ha) — 113 A. Such 
high equivalent width is only attained for stellar populations 
younger than ~6 Myr, which is consistent with the lower 
end of the age range obtained for Ap B from the continuum 
modelling. 

IRAS 17179+5444: a single 5 kpc aperture was ex- 
tracted for the analysis of the stellar populations in this 
single nucleus object. The optical spectrum of the nucleus 
is that of a Sy2 galaxy. The modelling results are consistent 
with a wide range of VYSPs ages, reddenings and percentage 
contributions. 

IRAS 20414-1651: this galaxy is classified as a sin- 
gle nucleus system in the new H S T H- band imaging sur- 
vey carried out by IVeilluex~et al.l ||2006|). Q pticallv classi- 
fied as an HII galaxv (|Kim fc Sanderslll998h . IRAS 20414- 
1651 shows no e vidence for a buried AGN at near- and mid - 
IR wavelengths (llmanishi et al. I. I2007I ). iFarrah et al] (|2003l ) 
modelled optical and IR photometric data and found that 
almost the entire IR emission in this source is powered by 
the starburst. Our modelling results reveal the presence of 
an important contribution to the flux from relatively highly 
reddened VYSPs. The age of ~6 Myr estimated using the 
Ha equivalent width (EW(Ha) — 43 A) is consistent with 
the lower end of the age range obtained from the modelling. 

IRAS 21208-0519: this double nucleus system is clas- 
sified as an Hll-galaxy in the optical. However, the northern 
nucleus, the brightest at optical and near-IR wavelengths, 
is the only region of the galaxy clearly showing Hll-like fea- 
tures in our spectra. The modelling results found for that 
nucleus (5kpc-I) are consistent the presence of moderately 
reddened, dominant VYSPs. The Hq equivalent width found 
for that region of the galaxy (EW(Ha) = 45 A) indicates 
VYSPs of age ~6.6 Myr. This is consistent with the lower 
range of ages found from the modelling for this component. 
However the modelling results found for southern nucleus 
(5kpc-II), reveal the presence of a dominant 12.5 Gyr OSP 
plus an "old" VYSP (90 - 100 Myr) with low reddening, 
that makes a small contribution to the fiux in the optical. 
Intriguingly, an important contribution from a VYSP com- 
ponent is again required in order to model the extracted 
spectrum from Ap B, centred in the bright blob to the south 
of the galaxy (see iKim et al.l . [20021 . their Figure 1). Finally, 
the modelling results for Ap A, sampling the region between 
the nuclei, are consistent with a mix of stellar populations 
from the two nuclei. 

IRAS 21219-1757: this object is classified as a Syl 
galaxy at optical wavelengths. Therefore, due to the power- 
ful AGN emission, it was not possible to model the stellar 
populations of this compact ULIRG. The extraction aper- 
ture and the extracted ID-spectrum are shown in Figures [2] 
and[6l 
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Figure 11. Extracted spectra from Ap B (upper panel) and Ap C 
(lower panel) for IRAS 22491-1808. The figure shows the detection 
of the He I absorption lines at 3820 A, 4026 A and 4388 A. These 
lines are not observed in stellar populations older than 100 Myr. 



IRAS 22491-1808: several knots and condensations 
located in both the tidal features and the nuclear region 
are visibl e in the high - resolu tion BIHK'-band images pre- 
sented in ISurace et al.1 (|2000l ). The estimated ages cover a 
wide range from 5-7 Myr up to a fe w hundred Myr. Wit h 
the optical spectrum of an HII galaxy (|Veilleux et al.l . [19951 ). 
IRAS 22491-1808 shows no signs of obscured nuclear activ- 
ity at IR wavelengths (Genzel et al., 1998; Veilleux et al., 
1999b; Imanishi et al., 2007) and it is therefore classified as a 
starb urst-powered ULIRG (|Lutz et al.l . Il999l : lArmus et al.l . 
|2003). Due to the small nuclear separation, and the fact 
that the slit is not perfectly aligned with the double nucleus 
structure, the 5 kpc aperture includes emission from both 
nuclei. Therefore, we decided to use another two apertures 
to sample the eastern (Ap B) and western (Ap C) nuclei 
individually. 

The modelling results found for both nuclei are consis- 
tent with the presence of moderatelly reddened VYSPs that 
make an important contribution, or even dominate, the fiux 
at optical wavelengths. Unsurprisingly, the modelling results 
for the 5 kpc aperture are consistent with the mix of stel- 
lar population in the two nuclei. The VYSP age determined 
from the Ha equivalent widths ( 64 A and 46 A for Ap B and 
C, respectively, corresponding to VYSPs of age ~6.5 Myr), 
is consistent with the results obtained from the modelling. 
Furthermore, the extracted spectra for Ap B and C (and 
also 5kpc) show the clear presence of He I absorption fea- 
tures at various wavelengths, which are indicated in Figure 
111! These lines are not observed in stellar populations older 
than 100 Myr and, therefore, provide further evidence for 
the presence of significant VYSPs in the two nuclei of the 
galaxy. Finally, Ap A was extracted to sample the tidal tail 
towards the north-west of the galaxy. Overall we do not find 
clear evidence for the presence of significant age/reddening 
gradients across the body of the galaxy. 

IRAS 23060+0505: this single nucleus system is clas- 
sified as a Sy2 in the optical, and show s clear signat u res of 
AGN activity at mid-IR wavelengths l|Genzel et al.l . Il998l : 
iLutz et al.l . ll999l ). For the nuclear region adequate fits are 
only obtained using Comb II, including a 10 Myr VYSP, a 
"red" power law component (p.l. index a — 6.1) and a 12.5 
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Gyr OSP making contributions to the optical light of 72%, 
6% and 22% respectively. It is worth remembering that these 
are percentage contributions in the normalising bin, which 
in the case of IRAS 23060+0505 is located at a wavelength 
range of 4450 - 4550 A. The contribution of the power-law 
component becomes more important with increasing wave- 
length, and is as high as 43% at 6800 A. This is consistent 
with the results of lVeilleux et al.l (|l997l ) who found a strong 
broad Pac line at near-IR wavelength, providing evidence 
for a significant contribution from a highly obscured AGN. 
For the other two apertures, Ap A and B, adequate fits are 
obtained using all combinations (i.e. the YSP properties are 
relatively unconstrained) . 

IRAS 23233+2817: iKim et all (|2002l l described this 
object as: "the only object in the 1 Jy sample with no ob- 
vious signs of interaction". Our Ap A samples a region to 
the so uth of the nucleus , including the bright knot detected 
m the iKim et al.l (j2002l ) K-band image. These authors de- 
scribed the latter as a bright HIT regions on the basis of their 
luminosities and colours. However we identified this knot as 
an M star in the galaxy based in our spectroscopic data (see 
Figure |6}. In the case of the 5 kpc aperture, adequate fits 
were obtained using Comb I and Comb 11, but not Comb 
III. Since the galaxy is classified as a Sy2, we concentrate 
on the results of Comb II. We find YSPs in combination 
with a 12.5 Gyr OSP, and power-law that makes a signifi- 
cant contribution to the flux at optical wavelengths. These 
results suggest that the AGN component is important in the 
case of IRAS 23233+2817. 

IRAS 23234+0946: this double nucleus system is 
classified as a LINER on the basis of its optical spec- 
trum. The western source is the brightest at optical wave- 
lengths, and shows no sign of AGN activity in the mid-IR 
ijlmanishi et al.l . l2007h . The modelling results found for this 
nucleus (5 kpc-I) are consistent with the presence of VYSPs 
that dominate the flux in the optical. Comparing these re- 
sults with those of the other two apertures, it is notable that 
the percentage contribution of the VYSP component in the 
western nucleus (60 - 82%) is significantly higher than in 
the eastern nucleus and the extended structure to the east 
of the galaxy (17 - 37%). Intriguingly, while high reddening 
values are found for the VYSP component {E{B — V)— 0.5 
- 1.5) to the east of the galaxy (5kpc-II and Ap A), lower 
reddening values {E{B — V) = 0.4 - 0.6) are found for the 
western source. 

IRAS 23327+2913: initially, two 5 kpc apertures 
were extracted sampling the two nuclei of this double nucleus 
system. In order to look for possible variations of the stellar 
populations within the galaxy when sampling larger regions, 
two additional, wider apertures were extracted, also centered 
on each nucleus. The extracted spectra for these apertures 
are identical to those of the 5 kpc apertures and therefore are 
not shown in Figure [S] In the case of the southern source 
(5kpc-I) an OSP is always required in order to model the 
extracted spectrum. The modelling results of Comb I are 
consistent with a mix of a 12.5 Gyr OSP plus a moderately 
reddened VYSP of 9 Myr that dominates the flux at optical 
wavelengths. On the other hand, this object is classified as 
a LINER at op tical wavelengths and shows no evidence for 
a buried AGN (jlmanishi et all . 120071 ). Therefore, it is likely 
that the power-law component used in Comb II accounts for 
a VYSP in this case, and it is possible that we are detecting 



the contribution of three different stellar components to the 
optical light: a 12.5 Gyr OSP plus important contribution 
from a highly reddened VYSP younger than ~ 5 Myr (rep- 
resented by a power-law with power-law index a = 1.54) 
and a VYSP of age 10 Myr and lower reddening. In the 
case of the northern galaxy (5kpc-II), it is possible to model 
the extracted spectrum using a single OSP of age 12.5 Gyr, 
with no need for any additional stellar component. Again, 
an OSP is always required to model the data. This demon- 
strates large differences between the stellar populations in 
the two nuclei. 

IRAS 23389+0300: the northern source of this dou- 
ble nucleus system is the brightest galaxy in the radio, 
known to hav e a milli-arcsecond radio core and extended ra- 
dio emission l|Nagar et al.l . l2003l ). Two 5 kpc apertures were 
extracted for this Sy2 galaxy, centred on the two nuclei. 
In order to perform a more detailed analysis of the stellar 
populations, an additional set of four smaller apertures was 
extracted sampling various regions of the galaxy. The ex- 
tracted spectra from Ap B and D are identical to those of 
5kpc-I and 5kpc-II, and hence, they are not shown in Fig- 
ure[6] The modelling results suggest the presence of younger 
VYSPs towards the north of the system, although we can- 
not rule out the idea of VYSPs of similar age at all locations 
sampled by the various apertures. In terms of reddening, we 
find no clear evidence for the presence of significant varia- 
tions across the body of the galaxy. 



